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ABSTRACT 


Probabilistic conposite micronechanics methods are developed that 
sinulate expected uncertainties in unidirectional fiber conposite 
properties. These methods are in the form of computational procedures 
using fbnte Carlo simulation. The variables in which uncertainties are 
accounted for include constituent and void volume ratios, constituent 
elastic properties and strengths, and fiber misalignment. A 
graphite/epoxy unidirectional conposite (ply) is studied to demonstrate 
fiber conposite material property variations induced by random changes 
expected at the material micro level. Regression results are presented 
to show the relative correlation between predictor and response 
variables in the study. These computational procedures nake possible a 
formal description of anticipated random processes at the intraply 
isuel, and the related effects of these on conposite properties. 
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CHAPTER I 


INTRODUCTION 


A. Background 

•Hie diverse requirements of recent engineering applications have 
motivated designers to explore specialized structural and naterial 
systems. Ceramic materials, for exanple, have several attractive 
structural properties, such as their high stif fness/weight ratios, and 
low variation of stiffness and strength over wide raixyes of 
environmental conditions. A significant disadvantage inherent to 
brittle structural materials is their vulnerability to failure die to 
cracks propagating from flaw. The increased probability of a flaw 
occurring in a naterial as the volume increases leads to bulk strengths 
idiich are a fraction of the theoretical strength of the naterial. The 
size effect on naterial strength (Ref. 1) can be explained by the 
"veakest link" concept. Griffith ( Ref. 2) reasoned that very snail 
solids, for exanple wires or fibers, might be expected to be stronger 
than large ones, due to the additional restriction on the size of the 
flaws. In the limit, a single line of nolecules mist possess the 
theoretical molecular tensile strength of a material. A consequence of 
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the size effect on strength ms the development of fiber conposite 
materials vfiich consist of thin v strong fibers bound together by a 
ductile natrix. The advantages of fine, strong fibers can explain the 
current trend tovaard increased use of fiber conposite materials in 
denanding aerospace applications. 

Properties of a conposite laminate depend on the properties of the 
constituent neterials, their distribution, and orientation. Laminates 
are composed of layers of unidirect ionally reinforced plies (laminae). 
The lamina is typically considered the basic unit of material in a 
conposite structural analysis, vtiich requires knowledge of the material 
properties of each individual lamina and its geonetric orientation. The 
branch of conposite mechanics that predicts ply material properties 
based on the properties, concentrat ion, and orientation of its 
constituents is known as conposite micromechanics, and frequently 
incorporates the traditional Mechanics of Materials assumptions. The 
desired laminate is created by stacking of plies in specific directions. 
The integration of ply properties to yield laminate properties is called 
laminate theory. Laminate variables such as ply orientation and 
stacking sequence can be tailored to yield a laminate with the desired 
material properties. Thus, the laminated conposite is a suitable 
material for conponent design. 

Analysis of fiber conposite structures is currently performed using 
a variety of computer codes. From the original codes based on classical 
micronechanics and laminate theory, recent codes (Ref. 3,4) have been 
developed tfiich incorporate the current state of the art. Complete 
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mechanical, thermal, and hygral properties are calculated, an! can be 
used to compute response. Advanced failure criteria are used to 
calculate composite strengths. Environmental effects are also 
quantified. The usefulness of these codes has been demonstrated by 


comparison with experimental and finite element results (Ref. 5,6). 

The analytical capability of many codes is limited by the 
deterministic nature of the computations. Specifically, fixed values 
for constituent material properties, fabrication process variables (i.e. 
const itient volume ratios) and internal geometry must be used as input. 
However, random variations in these parameters are not -only expected, 
but easily observed experimentally. (See Fig. 1) 

The analysis of composite structures requires reliable predictive 
models for material properties and strengths. However, the prediction 
efforts have been complicated by inherent scatter in experimental data. 
Since uncertainties in the constituent properties, fabrication 
variables, and internal geometry would lead to ircertainties in tte 
measured composite properties, the question arises: 


How ntch of the "statistical" scatter of experimentally observed 
composite properties can be explained by reasonable statistical 
distribution of input parameters in composite micro mechanics and 
laminate theory predictive models? 


The increasing use of probabilistic methods in structural mechanics has 
been shovn to provide a more realistic depiction of structural response 
due to load variations. (Ref. 7) The recognition that material 
parameters are characterised by a spectra of valves ( that is, are 
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statistical in nature 
probabilistic methods 


) rather than by a unique set of values, points to 
as a logical analysis approach. 
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B. Purpose 

The aim of this thesis is to develop a computational oapability to 
simulate the probabilistic variations in the mechanical behavior of 
unidirectional fiber composites. The Ibnte Carlo netted is used to 
simulate a variety of random processes, to quantify fiber conposite 
naterial variations induced by random changes in conposite fiber 
alignment, constituent properties, and fabrication process variables. 
This random process description is an attenpt to nore accurately predict 
the behavior of manufactured materials, which inherently include these 
random variations. The characterization of fiber reinforced conposites 
through sinulation of local nonuni formi ties provides an economical 
alternative to experimentation to measure material properties. 
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C. Formulation of the Mbdel 

The TTExiel oomonly used in characterizing fiber composites is based 
on the calculation of properties of the basic unit of an orthotropic 
ply. The layuqp geonetry is then used in laminate equations to calculate 
composite properites (See Figs. 2a, 2b). In this work, towever, the 
basic unit is taken as the sub-ply, vhich consists of only one 
fiber-matrix level in the material. Micro mechanics theory is used to 
calculate the properties of the assumed orthotropic sub-ply, each with 
randomly distributed fabrication variables and material properties. 
Distributed fiber directions, due to possible misalignment within the 
ply, are then used in the laminate equations to calculate ply 
properties. This sub structuring of the composite ply represents a novel 
attempt at characterization of fiber composite material properties based 
on probabilistically distributed constituent properties, individual 
fiber misalignment , and fabrication process variables (See Figs. 3a, 3b). 

This formulation is particularly well suited to the probabilistic 
description of fiber composite material properties. Since the 
micro mechanics and laminate equations can be used to calculate ply 
properties at any number of points in a ply, a tractable finite element 
structural analysis based only on simple distributional assumptions for 
physical parameter variations can be performed. This model supplies a 
rational procedure for composite material property assessment, because 
it treats the material as the result of a series of random processes 
which occur at the intraply level. 


(a) orthotropic ply 


(b) laminate 


Fig. 2- Conventional Model 



(a) subply 



(b) ply 


Fig. 3- Substructure Model 




9 


D. ffethod of Investigation 

!• ft'ief Description of ICAN 

•Hie Integrated Composite Analyzer (ICAN) is a computer program for 
comprehensive linear analysis of multilevel fiber composite structures. 
The program contains the essential features required to effectively 
design structural components made from fiber composites. It now 
represents the culmination of research conducted since the early 1970» S| 
at the National Aeronautics and Space Administration (NASA) Lewis 
Research Center (LeRC), to develop and code reliable composite mechanics 

theories. This user friendly, publicly available code incorporates 
theories for 


1 . 

2 . 

3. 

4. 

5. 

6 . 

7. 

8 . 
9. 

10 . 

11 . 

12 . 

13. 


conventional laminate analysis 
intraply and interply hybrid composites 
hygrai, thermal, mechanical properties anl response 
ply stress-strain infltence coefficients 
microstresses and mucrostress infltence coefficients 
stress concentration factors around a circular hole 
predictions of delamination locations around a circular tole 
oisson's ratio mismatch details near a straight free edge 
©age interlaminar stresses 
laminate failure stresses 


normal and transverse shear stresses 

explicit specification of matrix-rich interply layers 

finite element material cards for HASTRAN, KffiC 


A detailed description of ICAN can be found in Reference (3). The 
ICAN code and documentation are available throi^h COSMIC, the Computer 
Software Management and Information Center, Suite 112, Barrow Hall, 
Athens GA, 30602. 
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2. Sumrary of Variables 

The variables studied in this vork can be separated into tvo 
categories. The independent variables to be similated using random 

sampling consist of the following (see Fig. 4a for fiber coordinate 
system) j 


Geonetry: 

fiber orientation ai^le 
Fabrication variables: 
fiber volune ratio 
void TOlune ratio 
Fiber properties 

longitudinal elastic modulus 
transverse elastic nxxlulus 
shear modulus, 1-2 plane 
shear modulus, 2-3 pi are 
fiber tensile strer^th 
fiber conpressive strength 

ffetrix properties 

elastic modulus 
netrix tensile strer^th 
matrix conpressive strength 
matrix shear strength 


(THETA) 

(FVR) 

(WR) 

(EFP1) 

(EFP2) 

(d : P12) 

(GFP23) 

(SFPT) 

(SFPC) 


(E») 

(SWT) 

(SWC) 

(SMPS) 


The dependent variables to be calculated using I CAN consist of the 

following ply properties, measured about the neterial axes (see Fig. 
4b): 


nor no 1 modulus in 1-1 direction 
normal modulus in 2-2 direction 
shear modulus in 1-2 plane 
Poisson's ratio for strains in 2 
by stresses in 1 direction 
Poisson's ratio for strains in 1 
by stresses in 2 direction 
Coefficients of theriml expansion 
in 1-1 direction 
in 2-2 direction 
coupling coefficient 


(ECU) 

(EC22) 

(EC12) 

direction induced 
(NUC12) 

direction induced 
(NUC21) 


(CTE11 

(CTE22) 

(CIE12 



(a) fiber (b) material 

Fig. 4- Coordinate Systems «. 



Fig. 5- Order of ICAN input data cards 
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strengths in naterial directions 


longitudinal tensile (SCJOTT 
longitudinal conpressive (SCXXC 
transverse tensile (SCYYT 
transverse conpressive (SCYYC 
in-plane shear (SCXYS 


The descriptions above should be consulted periodically for the 
definitions of variables that henceforth will be referred to 
synbolically. 



3. Monte Carlo Methods 

Complicated stochastic processes can be s inula ted by a variety of 
numerical methods generally referred to as fbnte Carlo methods (Ref. 8). 
The term refers to that branch of experimental rat hematics concerned 
with experiments on random numbers. Since the advent of high speed 
computers , they have found extensile use in most fields of science and 
engineering, in analyzing many physical processes of a statistical 
nature, or where direct experinentation is not feasible. In general, 
they can be economically used to achiew a level of precision between 90 
and 95 percent. 

A Monte Carlo experinent refers to the procedure of randomly 
assigning a value to an independent random variable in a chosen model, 
and observing the dependent variable at the conclusion of the process 
being modeled . A Monte Carlo procedure is conposed of n such 
independent experiments, then n is sufficiently large, the observations 
will yield, by virtue of the lavs of large numbers, a statistically 
meaningful description of the physical problem. 

The form of Monte Carlo used in this study is as follows: 

1. Define the system model by assuming 

a. model regression function 

b. method of error incorporation 

c. probability distributions of all errors (for all independent 
variables) 

d. any eqiations used to model the phenomena of interest 

2. Use the computer and random sampling techniques to select 
values of the independent variables. 

3. Calculate dependent (output) variables using the prescribed 
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equations. 

4. Est irate regression parameters for the assumed model. 

5. Replicate the experiment, each tira with a iew set of input 
values. 

6 Use appropriate statistical nethods to calculate properties of 
the distribution of parameter estirates. 


E. Brief Summary of Results 

A ply nade from the AS-Graphite /IMS epoxy composite system is 
studied. Tte monte carlo scheme is used to generate a lumber of 
response results, vhich are analyzed in graphical and numerical form, to 
supply a random process description of composite ply elastic constants, 
t terms 1 expansion coefficients, and strengths. Histogram and 
distribution plots of results for assumed narrow and wide variations in 
input properties are compared with a deterministic base case for an 
aligned ply. The figures demonstrate the range of values that response 
variables assume for the exanple data under consideration. 

Confidence intervals are calculated for response variables in 
subsequent samples, which are normalized with respect to an appropriate 
independent variable, to yield plots of normalized response as a 
function of fiber volume ratio, for various valiss of distribution 
parameters for the related independent variable. These plots 
demonstrate tte sensitivity of ply properties to randomly selected 
uncertainties in constituent and fabrication variables. 

Several multiple lirraar regression models wre calculated for 
response variables. Tte relative correlation of predictor (independent) 
variables with response is studied for all output properties considered. 
Varying levels of significance were achie**d in tte regression 
equations, due to tte differences in complexity of response variables. 
Elastic constants can be described adequately with simple regressor 
functions, and generally explain betveen 80 and 99 percent of tte 
observed response variations about a mean. Tte regression models 



16 


studied for strength, although achieving better reliability with higher 
order regressor functions, demonstrate such low signifioanoe as to be 
practically useless for predictive purposes. This is rot an unexpected 

re5U ^i because of the conplex nature of strength behavior in composite 
materials. 



CHAPTER II 


rETHODG OF CALCULATION 


A. Overall plan 

1. Input stricture for I CAN 

The input data for a typical execution of the available I CAN 

program consists of (see Fig 5) 

1 . header card 

2. control cards 

3. ply data cards 

4. naterial system cards 

5. load cards 

For repeated use of the I CAN program, input data files nust be 
created and used one at a tine. Each successive run of the master 
program (of which I CAN is Trade a subroutine) vrites the input file from 
user-supplied paraneters and calls I CAN. The ply data cards contain 
randomly generated fiber orientation angle values. The material system 
cards contain randomly generated values for fiber and void volume 
rat ios . 
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2. Constituent Property Variations 

Each successive execution of ICAH uses a distinct set of material 
properties for fiber and matrix. Hie random nunher generation is 
perforned with user-supplied parameters which are stored in a separate 
file. The options of using either generated properties or using the 

contained in the resident data bank are available. Any subset of 
the parameters described nay be generated or held constant with proper 
specification of the Boo leans which control the input to the I CAN 
program, (see Figs. 6,7) 
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^555*57®^®^ VARIES; CONSTANT FIBER VOLUME RATIO 


STDATA 15 1 

T 50 T F 
F 000.0 10.0 

F 
F 
T 

PLY 

MATCRDAS-1IMHS 
PLOAD 10. 0.0 

PLOAD 0.0 0.0 

PLOAO 0.0 0.0 

OPTION 0 


15 
T 

0.300 


T 
T 

0.200 


70.00 70.00 

AS-1IMHS 
0.0 0.0 
0.0 


3.00 


.0 

0.0 


OF 0 . 30 ; TAPE 003131 
5 


.000 

.57 .03 


Fig. 6- Command Input 


EFPl 

T 

0.3100E 

oa 

0.3000E 

07 

EFP2 

T 

0.2000E 

07 

0 . 2000E 

06 

GFP12 

T 

0.2000E 

07 

0.2000E 

06 

GFP23 

F 

0.1000E 

07 

0. 1000E 

06 

SFPT 

T 

0.4000E 

06 

0. 1000E 

02 

SFPC 

T 

0.4000E 

06 

0.1000E 

02 

EMP 

T 

0.5000E 

06 

0 . 5000E 

05 

SMPT 

T 

0.1500E 

05 

0.1000E 

02 

SMPC 

T 

0.3500E 

05 

0.1000E 

02 

SMPS 

F 

0.1300E 

05 

0.1000E 

02 


Fig. 7- Constituent Variation Input. Example for AS-1 Graphite 
fiber and IMHS Epoxy matrix, with wide variations of 
stiffnesses and strengths. 


20 


3. Repeated rims 

The user nust specify the number of I CAN runs desired in a given 
sample. In this study, fifty (50) runs were used throughout, to take 
advantage of the simpl if ication in statistics by using suitably large 
sanples. From elementary statistics, it is known that any process that 
is the result of the combined interaction of several probabilities can 
assumed to approximate a normal distribution. For phenomena that are 
assumed to approximate a normal distribution, the simplest forms for 
calculating statistics apply to suitably large sanples (usually greater 
than thirty). The sample size of fifty vas chosen to supply a 
practicably large amount of data, within the restrictions imposed on 
conputat ion tine. 

The data generated by repeated execution of the ICAN routines is 
stored in a sequential access dataset, where the 50 output files are 
separated by end of file markers. This arrangement allows a single 
Fortran unit to be used for output throughout. A simple flowchart of 
the data generation routines is shown in Fig. 3(a). 

4. Data collection 

The ICAN output files are searched to locate the specific material 
properties and strengths of interest in this study. The flowchart of 
data collection routines is shove in Fig. 8(b). After obtaining the 
sample of ICAN output, the investigator may choose to scrutinize 
parameters or calculate statistics aside from those chosen in this 
study. This is likely, in light of the large quantity of data available 
and the need for limiting the scope of this particular study to 
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representat iue pro p erties. The user would haw to supply additional 
code or adapt existing code to suit his purposes in this case. The 
coded nod ificat ions to ICAN used in this study are included in Appendix 
A. 
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B. Generation of Pseudo Random Nuiribers 

An integral part of any ironte carlo simulation is the use of random 
nuiribers having a specified distribution vftich is assuned to characterize 
the process under study. Indeed, many statistics textbooks carry tables 
of raixlom nuiribers as appendices. Simulations using large sanples 
require many repeated calculations, each with different ''random” 
nunfcers. Sire* filling of a computer memory with a large table of 
random nuiribers is wasteful, algorithms have been developed (Ref. 9) to 
generate streams of random nunfcers riienewr needed in the process of 
calculations. The nuiribers used are usually obtained using some form of 
a recursion relation, hence the sequence is termed pseudo-random. 

1. Uniform Distribution 

The starting point for many random nuntter schemes is the uniform 

random nunfcer generator, rfiich simulates a sample from the uniform 

distribution. A continuous random variable has a uniform distribution 

over an interval a to b ( b > a ) if it is equally likely to take on any 

value in this interval. The probability density function is thus 

constant over ( a , b ) and has the form 

1 

f («) = “b^T a £ x - b 

= 0 elsevfiere 

The probability distribution function is, on integrating 

F(x) =0 x ( a 

x - a 


b - a 


a < x < b 
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= 1 x ) b 

The uniform distribution is shovn in density and distribution form in 
Figs. 9a and 9b. 

Le hirer (Ref. 1©) proposed the congruent ial method of generatirxj 
pseudo random nunbers conforming to the uniform distribution. Tte 
recurrence relation takes the form: 

x^ = (ax^_j + b ) modulo m 

vfiere the notation signifies that x. is the remairaler when (ax. , + b) 
is divided by m. The nultiplier a, increnent b, and nodulus m are 
integers. The starting value x^ nust be assumed, and is known as the 
"seed" of the generator. Generators for viiich b = 0 are krxsvn as 
nultiplicat iw. They are called mixed vben b is nonzero. Because 
selection of the nultiplier a and msdulus m strongly influence the 
generator, most generators in use are of the nultiplicatiue form. A 
discussion of the choice of parameters, naximum period, and degree of 
correlation of this generator is available (Ref. 11). 

For a given uniform random nunber u on the interval ( 0 , 1 ) a random 
nunber x having a desired distribution F(x) is often obtained by solving 
the equation u = F(x) for x (Ref. 12). Since the process requires the 
determination of the inverse distribution function F *(x), its use 
depends on the ease of deriving the expression or some approximation. 

The following sections describe the distributions used, and net hods for 
generating random nunbers on those distributions. 
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2. Normal (Gaussian) Distribution 

Tbe nost ooimon distribution is the familiar nomal distribution, 
with the "bell shaped*' density function, given by 



< k < «>, f* < «, and a > 0 

with mean ft and standard deviation a. The distribution function is 
written 


F(x) 




wfiich cannot be expressed in closed form analytically but can be 
numerically evaluated at any value of x. 

Tbe Bdx-MjI ler or "Polar" method (Ref. 13) is most conronly used 
for generating random deviates from a mean to approximate tte nornel 

distribution. If x^ and x^ are independent uniform ranlom variables, 
then 

y A = o(-2 In x^)®* 5 cos 2nx 2 ♦ ft 
y 2 ~ ff (“2 In x,)®* 5 sin 2 nx^ ♦ ft 


are independent random variables with the standard mrnal distribution 
having mean ft and standard deviation a. 
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3. Ganna Distribution 


He gairam distribution is a two-paraneter distribution whioh is 


flexible in fitting a variety of random processes. It is a one sided 
distribution in that physical quantities that are limited to values in 
the positive range are frequently modeled by it. Its density function 
is giwn by 


f(x) = 


n k) 


-Xx k-1 

S X 


where x, X, k > 0, and k is an integer. 

Tte paraneters X and k nay be interpreted as scale and shape parameters, 

respectively. J*(k) is the well known gantra function, 

. fco k-1 -u. 

^(k) = J 0 u e du » 

which is widely tabulated. The gamra distribution function is given by 


F(x) = 


n k) 


0 


k-1 -Xu . 
u e du 


r(k,Xx) 


x > 0 


' F(k) 
s 0 else viie re 


tftere r(k,u) is tte inconplete gantm function 

u 

r(k,u) = 

0 


k-1 -x. 

x e dx 


which is also widely tabulated. For integer values of k, 


T(k) = (k-l)» 

ani tte gamma distribution is know as the EWangian distribution after 
A. K. Erlang y vho introduced it in the theory of qietes and Markov 
processes. 
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Ganna variates are generated using the sequence 

U 1 ,U 2 ,U 3’ “k 

satisfying the uniform distribution on the interval (0,1). 
The recursion relation is 


y i " - — 


in u. , 

i* 


£ lk 

= ) y. = — in B u. 

A 1 x i=i 1 


1=1 

viiere x is a gamma variate having parameters X and k (Ref. 14). 
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4. Ueibull Distribution 

The (bibull distribution (Ref. IS) is nost popular viien modeling 
problems of reliability, imterial strength, and fatigue. The Ueibull 
density function is given by 

f(x;<x,£) = a/Jx^~*exp(-ax^) 

0 < x < eo, a > 0, ft > 1 

viiere a and ft are the shape and scale parameters, respectively. The 
cunulatiw distribution function 

y = F(x) = 1 - exp[-(x//J) a ] 

leads imnediately to the inverse relationship 

F _1 (y) = x = - ln(l-y) ] l/a 

as the desired Ueibull random generator *iien y is a uniform random 
variable. 

Figures 9-12 show the above distributions in analytical form. 
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C. Distribution Assumptions 

The variables chosen for variation are those for «hich reasonable 
assumptions can be rode to describe their distribution. The fiber 
geometric configuration with respect to ply axes is assured to follow a 
normal distribution with mean of zero (degrees) and sore srell standard 
deviation, to be specified. The fiber volume ratio is assumed to be 
norrolly distributed about sore mean between 0.3 and 0.7. The void 
volure ratio, which is ideally srell, is assured to follow a ganra 
distribution skewed toward zero. (Hote that in the gamma distribution 
used, a value of zero has a probability of zero. This model is chosen 
because the state of rest present renufacturing technology precludes the 
fabrication of a fiber conposite conpletely free of wid.) 

The properties of individual fibers and ratrix are varied. The 
normal and shear noduli are assured to follow the norrel distribution, 
and tte strengths are assured to be Ifeibull distributed. 

Figs. 13-27 show the results of random number generation in each 
distribution studied. The density (or histogram) and cumulative 
distribution plots are shown. Several veibull and gamma distribution 
simulations are shovn, to demonstrate the effects of assured parameter 
variations on the distribution sampling. 
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. HISTOGRAM FOR 

A HE I BULL GENERATOR 



.DISTRIBUTION OF 

HE 1 BULL GENERATOR 






HISTOGRAM FOR 
HE I BULL GENERATOR 


RANGE t ( Ot) (ksi ) 
(a) histogram 



6 « 400 ksi 
a = 15 


DISTRIBUTION OF 
HE 1 BULL GENERATOR 

B * 400 ksi. 


RAN8C tt on (ksi) 

(b) cumulative distribution 
Fig. 26- Weibull Distribution Simulation 

Fiber Tensile and Compressive Strength 
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HISTOGRAM FOR 
HE I BULL GENERATOR 



DISTRIBUTION OF 
HE 1 BULL GENERATOR 
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D. Use of I CAN 

•nils section describes the essential theories and assumptions 
incorporated in the ICAN program. The symbolic notation conventions, 
formulations, and definitions are incited in Appendix B. 


1. Composite Micro mechanics 

The branch of composite mechanics viiich relates ply properties to 
constituent properties is known as composite micronechanics . The inputs 
consist not only of constituent material properties (fiber and matrix), 
but geometric configuration and fabrication process. Output includes 
ply hygral, ttermal, and mechanical properties. The assumptions for 
equation development are: (Ref. 16) 

1. The Mechanics of Materials are used to derive the equations, 
allowing each property to be individually identified. 

2. The ply resists in-plane loads according to the schematic 
shown in Fig. 4(b). 

3. Tte ply and its constituents behave in a linear elastic manner 
to fracture (see Fig. 28). 

4. Tte ply is transversely isotropic in the 2-3 plane. 

5. Tie matrix is isotropic. 

6. Complete bond exists at the fiber-matrix interface. 

The direction conventions and terminology used in the equations 


are: 

1. Properties measured along fiber direction are called 
longitudinal. 

2. Properties measured transverse to fiber direction are called 

trams verse. . 

3 In-plane shear is also knovsn as intralaminar shear. 

4*. All ply properties are defined with respect to ply material 
axes (1,2,3) for description and analysis. 
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2. Laminate Theory 

Classical laminate theory supplies a oonwnient procedure to 
predict the response oF a laminate to external load. The theory uses 
anisotropic elasticity to obtain the stress-strain relationship For the 
basic lamina. The stress-strain relations oF individual laminae are 
transForned to coincide with a global set oF reFerence axes. The 
stress-stain law oF the laminate in terms oF the properties and 
distribution oF individual laminae are calculated using a sunration. 
Resultant Forces and moments are deFined by integrating the stresses 
through the thickness oF the laminate. The plate const itutive equation 
is inverted, giving midplane strains and plate curvatures in terms oF 
applied Forces and monents. These strains and curvatures are 
substituted into the lamina stress-strain equation to obtain lamina 
stresses in the global system. The stresses obtained are then 
transForned into the principal neterial system oF the lamina in question 
and conpared with ultimate stresses obtained using Failure criteria. 
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3. Strength Theories 

The strength theories in ICAN make use of sewral assumptions. 
First, it is assumed that there are five characteristic values of 
strength of a unidirectional conposite: 

1. longitudinal tensile strength 

2. longitudinal conpressive strength (3 separate criteria) 

a. rule of mixtures 

b. fiber microbuckling 

c. de lamination 

3. transverse tensile strength 

4. transverse conpressive strength 

5. in-plane or intralaminar shear strength 

He fracture nodes usually associated with these strengths are shown 
schematically in Fig. 29. 

Once ply strengths are calculated (in the ply coordinate systems), 
geonetric transf ornat ions are used to calculate conposite failure loads. 
The process used is briefly described below. 

1. Calculate loads (in conposite system) required to induce load equal 
to ply strengths (in ply systems) for each node. 

2. Calculate minimum of failure loads for each ply. 

3. Calculate minimum of failure loads of all plies, and use this load 
as the failure strength of the conposite for a particular failure 


node. 
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E. Review of Applicable Statistical Concepts 

Conposite properties are calculated for large sanples using a 
specific set of distributions of input properties. In this context, 
snail sanpling theory does not apply, because the sanples used are 
sufficiently large. 

1. Sanple Means 

Calculation of tte naan sanple values proceeds by defining 

n 

L x 

_ 1=1 1 

mean = x = — 

n 

vfiere n = sanple size 

x .= sanple data valves 

The population mean is unknown, so the sanple mean is assumed to be the 
best estimator of the population mean. 

2. Sanple Standard Deviation 

An estimate of the population standard deviation is calculated 
using tte statistically efficient estimator 

T n r - 1 1/2 

° m ("i * “i* j • “ 2 ” 

3. Confidence Interval Estimates 

An inportant problem in the area of statistical inference is the 
estimation of population paraneters (such as mean, variance, etc.) from 

Parameters x and a are the mean and standard 


sanple statistics. 
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deviation of tte sampling distribution of a statistic S. Hie sampling 

distribution of S is assumed as approximately normal (idiioh is true for 

nany statistical distributions if n > 30). Confidence interval 

estimates are constructed for the statistic S. Thus, intervals are 

identified for tfiich it can be asserted with a reasonable degree of 

certainty that they contain the parameter considered. Obviously, the 

degree of certainty (or confidence level) will vary with the size of the 

interval chosen. Values of confidence coefficients, z^ , are associated 

with confidence levels. For example, an actual sample statistic S is 

expected to be found lying in the interval (x - z c a) to (x + z c o) (where 

a is tte unkreiwi population standard deviation) some percent of the 

time. Let the z value in this exanple be 1. Assuming a normal 
c 

sampling distribution, (with z^ = 1) the normal distribution area 

function specifies that S falls between (x - a) and (x + o) about 

68.27/. of the time. Similarly, the confidence of x lying in the 

interval (S - a) to (S + a) is about 68.27/.. The endpoints of the 

intervals are known as confidence limits. Various confidence 

coefficients z , corresponding to frequently used confidence levels, 
c 

have been tabulated. 

In this work, the confidence interval for means is given in terms 
of the sample statistics by 

a 

x + z — 

" c ST 

vtere z is the confidence coefficient, vfiich takes on values of 
c 

1.645, 1.960, and 2.580 for the 90, 95, and 9971 confidence levels, 
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respect ively. 

4. Regression 

The term "regression" as used in the area of statistics refers to 
the process of formilating a mtheratical model to explain randomly 
observed phenomena. So ite functional form for the vay each variable 
enters the model must be assumed. Comparison of the degree of fit of 
different assumed nodels ideally leads to a better model. The basic 
regression strategy used here consists of : 

1. Assume a multiple linear regression model. The normal equations 

for such a model are: 

{V) - [X]<0> + {0 

viiere 

{Y} = vector of dependent variable values 
[X] = matrix of functions of independent variable 
= regression "true" values 
= errors 

The normal equations can be solved as follows: 

[X] T {Y> = [X] T [X]{/J> + [X] T {e> 

{b> = [x t x] _1 [xT{y> 

viiere 

{b} as parameter estimates 

2. Use a standard statistical package (Ref. 17) to estimate regression 

parameters- 
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3. Calculate properties of regression parameter distributions to 
assess model precision. 

T 

In the event that [X X] is singular, inplying that some of the 

T -1 

normal equations are linearly dependent, [X X] does not exist. The 
model should be expressed in terms of fewer parameters, or stould 
include assumed restrictions on the parameters. 

2 

The square of the multiple correlation coefficient, R , is usually 
calculated for each regression model, and supplies a convenient measure 

A 

of the degree of fit between data values {Y} and values {Y} 
predicted by the regression equation. It is defined by 


R 


2 


Sum of Squares due to regression model 
Total Sum of squares about mean Y 


L (Y. - Y) 2 
L (Y. - Y) 2 


Frequently, it is necessary to determine if inclusion of particular 
terms in a regression nodel is worthwhile. To this end, the extra 
portion of the regression sum of squares which arises due to the terms 
under consideration is calculated. The mean square (defined as the sum 
of squares divided by the corresponding degrees of freedom) derived from 
this extra sum of squares can be conpared with s 2 , the estinate of o 2 , 
to see if it appears significantly large. If it does, the terms under 
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consideration should be included. The statistic is frequently compared 
to the appropriate percentage point of the F- distribution, which is 
tabulated. 

Supopose the extra sum of squares due to a paraneter, given that a 
number of other paraneters are already in the model, is calculated. 
Symbolically, 


SS ( b il b 0' b l»* , -» b i_l* b i+1 i***»l» k ) i = 1,2, ...,k 
represents a one degree of freedom ( 1 df ) sum of squares viiich 


Treasures the portion of the regression sum of squares due to the 
coefficient b^ This is a measure of the value of addirp a term to 
the model viiich previously did not include 0.. The corresponding mean 
square, equal to the SS (since it has one df) can be compared by an 
F- test to s 2 . This is knovn as a partial F- test for the single 
parameter viiich is a special case of the F- test described earlier. 

The stepwise regression procedure (Ref. 18) is a structured way to 
insert variables in order of correlation until the regression equation 
is satisfactory. The partial correlation coefficient measures the 
relative importance of terms not yet in the model, to choose the next 
candidate for entry. The analagous statistic, F- to enter (or F- to 
remove) is usually evaluated for each predictor at every stage as though 
it vere the last term to enter the model, to determine if terms retained 


at a previous step have become superfluous, because of some linear 
dependence with terms now in the model. The largest F— statistic 


calculated at each step is compared with the appropriate percentage 
point of the F- distribution, and the predictor variable is entered (or 
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relieved) based on the significance of this F- test. Testily of tte 
least useful predictor is perforned at every step. The R 2 statistic is 
calculated, to provide a Treasure of the value of the regression at each 
step. This stepwise linear regression scheme is used in this work 
because of its computational economy, and because it allows the analyst 
to assess the relative influence (or correlation) between individual 
predictor variables of a selected model and response for a particular 
data sample. Other schemes are available (Ref. 18), such as backward 
elimination. The stepwise procedure is recommended for its direct 
nature in testing the model with only significant predictor terms. 



CHAPTER III 


RESULTS 

ft- Property Histograms and Distributions 

In this vork, fiber and matrix properties are allowed to assune a 
range of values to assess the sensitivity of the conposite ply 
properties to constituent perturbations. Graphite fiber and epoxy 
natrix are used as the constituents. Initially, two separate samples of 
output data are generated and studied to demonstrate the effects of 
input parameter changes on conposite material properties. Ttese two 
cases are conpared with a deterministic base case with no random input 
property generation. The data for all three cases is given in Table I. 

The results of cases 2 and 3 are shown in histogram and cunulative 
distribution form in Figs. 3® - 42. The results of tte deterministic 
case 1 are summarized in Table II, and can be easily conpared with the 
histograms and distributions. 
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TAH-E I- 

INPUT DATA FOR SAMPLING 


INPUT 


CASE 1 

CASE 2 

CASE3 

THETft ( degrees ) 

0.0 


_ 

A* 


- 

®.® 

0.0 

a 


- 

5.0 

10.0 

FVR 


0.50 

- 


A* 


- 

0.5 

0.5 

a 


- 

0.1 

0.2 

WR 


0.01 

- 


X 


- 

3.0 

3.0 

k 


- 

3 

5 

EFPl(ksi) 


31000 

- 


A* 


- 

31000 

31000 

a 


- 

1500 

3000 

EFP2(ksi ) 


2000 

- 


V 


- 

2000 

• 2000 

a 


- 

100 

200 

GFP12(ksi ) 


2000 

- 


A* 


- 

2000 

2000 

a 


- 

100 

200 

Q r P23(ksi ) 


1000 

- 


t* 


- 

1000 

1000 

a 


- 

50 

100 

SFPT(ksi) 


400 

- 


£ 


- 

400 

400 

a 


- 

20 

10 

SFPC(ksi) 


400 

- 


0 


- 

400 

400 

a 


- 

20 

10 

EJF(ksi) 


500 

- 


A* 


- 

500 

500 

a 


- 

25 

50 

SIf>T(ksi) 


15 

- 


fi 


- 

15 

15 

a 


- 

20 

10 

SfFCfksi ) 


35 

- 


fi 


- 

35 

35 

a 


- 

20 

10 

SPPS(ksi) 


13 

- 


fi 


- 

13 

13 

a 


- 

20 

10 


TABLE II- CASE 1 RESULTS 


PROPERTY 


VALUE 


ECU 

EC22 

EC12 

NUC12 

NUC21 

CTE11 

CTE22 

CTE12 

SCKXT 

SCXXC 

SCYYT 

SCYYC 

SCXYS 


15750 ksi 
1065 ksi 
516 ksi 
0.275 
0.018 

0.775 x 10 H 
0.181 x 1 0 _< 
0.000 
203 ksi 
165 ksi 
11.74 ksi 
27.41 ksi 
10.01 ksi 
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HISTOGRAM FOR ECU HISTOGRAM FOR ECU 



DISTRIBUTION OF EC11 DISTRIBUTION OF EC1 1 



CUMULATIVE FREQUENCY 








CUMULATIVE rBEOUEMCT 



(c) case 2 distribution (d) case 3 dls 

Fig. 32- Sampling results for In-plane Shear Modulus 




CUMULATIVE FREQUENCY 








CUMULATIVE 
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HISTOGRAM FOR CTEll 

long, thermal expansion 


HISTOGRAM FOR CTE 
LONG. THERMAL EXPANSION 
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HISTOGRAM FOR CT 122 


HISTOGRAM FOR CTE22 



2 


DISTRIBUTION OF CTE22 DISTRIBUTION OF CTE22 



Fig. 36- Sampling results for Transverse Thermal Expansion 



CUMULATIVE 






FREQUENCY 


histogram for scxxt 

10N6. TENSILE STRENGTH 



RANGE < E 03> 

(a) case 2 histogram 


2?. 5 

tei 2 • 0 


histogram for scxxt 

LONG. TENSILE STRENGTH 



• I e .16 ,20 

RANGE < E 03) 

(b) case 3 histogram 


distribution of scxxt 

LONG. TENSILE STREN6TH 


distribution of scxxt 

LONG. TENSILE STRENGTH 



RANGE ( E 03) 

(c) case 2 distribution 


RANGE ( E 03) 

(d) case 3 distribution 


Fig. 38- Sampling results for Longitudinal Tensile Strength 
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HISTOGRAM FOR SCXXC 


HISTOGRAM FOR SCXXC 



DISTRIBUTION OF SCXXC 


DISTRIBUTION OF SCXXC 



Fig. 39- Sampling results for Longitudinal Compressive Strength 
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HISTOGRAM FOR SCYYT 

TRAN. TENSILE STRENGTH 


HISTOGRAM FOR SC' 

TRAN. TENSILE STRENGTH 




mm 

f'VVimtiviaa 

kange < E 0?) 

(a) case 2 histogram 

RANGE ( 

(b) case 3 

E 02) 
histogram 

DISTRIBUTION OF SCYYT 

DISTRIBUTION 

OF SCYYT 

TRAN. TENSILE STRENGTH 

TRAN. TENSILE 

strength 


RANGE <E 02) RAN6E < £ 02) 

(c) case 2 distribution (d) case 3 distrlbuti 

Fig. 40- Sampling results for Transverse Tensile Strength 
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CUMULATIVE 
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HISTOGRAM FOR SCXYS 


HISTOGRAM FOR SCXYS 



DISTRIBUTION OF SCXYS 


DISTRIBUTION OF SCXYS 



Ffg. 42- Sampling results for In-plane Shear Strength 
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B. Fiber Strength Effect 

To show the effect of fiber strength change, on the longitudinal 
strength, of the composite, several shape parameter, of the eeibull 
distribution for fiber strength are assured. The rente carlo procedure 
is then condoled at several fiber Vo lure ratio values, fill properties 
are oaried, encept fiber volune ratio. The distribution paraneters of 
all properties eacept fiber strengths are held constant. The curves 
generated are shows in Figs. 43 and 44. In the figures the solid lines 
and syrhols show the mans of the 957. confidence internal estiretes for 
the sanple size of 50 chosen at each point. The points on both sides of 
each curve locate the upper and loser bounds of the confidence 
intervals. The convention described is intended to provide a convenient 
indication of the dispersion of the sanple valms at each point. 


NORMAL I ZEO (SCXXT/SFPT) 


LuNGo TENSILE STRENGTH 
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long, compress, strength 



Fig. 44- Longitudinal Compressive Strength; for various 
shape parameters of fiber strength. 
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C. Ifetrix Strength Effect 

The effects of changes in matrix strength on composite strengths 
are studied by suitable variation of the shape p^aseters governing the 
netrix strength distributions. Analagous to the plots given for fiber 
strength effects, the matrix effects are stovn in Figs. 45 - 47. 
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TRANSo tensile strength 



F*g. 45- Transverse Tensile Strengthi for various 
shape parameters of matrix strengths. 


NORMALIZED (SCYYC/SMPC) ( E -02) 
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TRANS. COMPR. STRENGTH 



Fig. 46- Transverse Compressive Strength; for various 
shape parameters of matrix strengths. 
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D. Fiber Orientation Effect 

“”'™ d V * 1 "” ° E th » orient.,, on distribution 

paraeeter are consecutively used in the nonte carlo procedure to assess 

the effects on several conposite properties. These plots are show, in 
Figs. 48-57. 


E. Fiber Stiffness Effect 

Assured ual«s of the fiber nodulus distribution parameter are used 
in the sinulation to similarly assess the effects on the related 
conposite properties. The plots thus generated are shove, in Figs. 


58-67. 


NORMAL I ZED (ECtl/ErPl) 



NORNAIIZEO (EC22/EMP) 


i MANS* ELASTIC MODULUS 


? 32 

tU _ 

v30 


A o « io e 

□ o » 5° 
V O * 1° 


1 ° Q ~ 7i 1 

* * • ' *5 ,6 

FIBER VOLUME RATIO 

Fig. 49- Transverse Elastic Modulus; for various 
shape parameters of fiber orientation. 
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FIBER VOLUME RATIO 
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NORMALIZED (ECI2/GFP12) 








NUC2I ( E _ 03) 
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POISSON’S RATIO (MINOR) 




NORMALIZED (SCXXT/SFPT) < E -02) 
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LONG* TENSILE STRENGTH 



Fig. 63- Longitudinal Tensile Strength; for various 
shape parameters of fiber modulus. 







NORMALIZED (SCYYC/SMPC) 




NORMALIZED <SCXYS/SMPS> ( E -01) 
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IN PLANE SHEAR STRENGTH 



Fig. 67- In Plane Shear Strength; for various 
shape parameters of fiber modulus. 
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G. Regression Models 

Hie output data of cases 2 through 11 are used as successive inputs 
to the regression scheme. The goal of stepwise regression, as used 
here, is to measure the degree of correlation between a dependent and a 
set of independent variables for a given set of data. Tte outputs of 
the regressions conducted show the independent variables accepted into 
the model (based on F-test criteria) in order of degree of correlation 
with the dependent variable of interest, along with the final R 2 
statistic. (The R valies represent the square of tie .multiple 
correlation coefficient, a convenient measure of the fit between data 
values and values predicted by the regression equation.) 

The ordering of predictor variables by stepwise regression has 
several important uses. In this study, the scheme facilitates easy 
investigation of the effects of material changes on composite 
properties. Since the monte carlo scheme permits generation of large 
amounts of data, the regression is easy, inexpensive, an! can provide 
insight concerning the sensitivity of dependent variables for assumed 
distributions of predictor variables. A variety of material 
configurations and constituent distributions are examined, and a model 
constructed for each dependent (or response) variable. It must be noted 
that the relative correlations of predictor variables with response 
variables will be functions of the assumed distributions, the particular 
data sample considered, and the functional manner in which the predictor 
variables are incorporated into the model. 

A simple regression model was assumed for each 


response variable. 
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The first set of "simple” regression models uses as predictor functions 
only the independent variables as individual terms. To be more precise, 
the predictor variables used are not simply the independent variable 
values, for there are 15 of these for each layup. The arithmetic mean 
of independent variable valws is thus used as the predictor variable in 
the first set of regression models. The only exception to this is the 
use of the sin 2 of the average of the fiber orientation aisles as the 
angular dependence predictor, denoted by THETA in the tables to follow. 
The simpler response variables can be adequately described using the 
linear function forms in the regression models. The simple variables 
include the elastic constants, (EJCll, E3C22, EJC12, NUC12, NUC21) and 
coefficients of thermal expansion (CTE11, CTE22) . The results of the 
regressions performed in the “simple" manner are given in Tables III - 
XIV. In the tables the input labeled with N1 through N5 and W1 through 
U5 represent narrow and wide distributions of all properties. Input 
labeled N6 through N10 and U6 through U10 describe the same 
distributions, except that the composite is assumed unidirectional, i.e. 
no angular variation. The distinction shows the reduction in predictive 
capability induced by deviations of the fibers from aligned orientation. 

The models assumed for the response (output) variables are of the 

Form 

Y = B 0 + BiX, + B 2 X 2 + B 3 X 3 + ... + B X 

n n 


where 


Y = response variable (ECU, EC22, E3C12, etc.) 
B = regression parameters to be obtaired 
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x„ - .«raae of independent nariable ual« s thresh the 
thickness of the ply (THETA, FVR, WR, . to ., 

^h " Dd " 1 POS, " 1 **' d “ -U independent variables that 

appear in the equations for the relate** n i 

he related ply property (see Appendix Bj . 
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TABLE III- LONGITUDINAL MODULUS fEClll 


SlfPLE fEDEL 


INPUT 

FVR 

TERTE ACCEPTED 

R a 

N1 

0.3 

FVR.EFPl 

83.17 

N2 

0.4 

FVR,EFP1 , THETA 

92.63 

N3 

0.5 

FVR.EFPl, THETA 

94.02 

N4 

0.6 

FVR, EFP1 , THETA 

94.59 

N5 

0.7 

FVR.EFPl 

84.00 

W1 

0.3 

FVR, THETA, EFP1 

64.49 

U2 

0.4 

FVR.EFPl , THETA 

89.88 

W3 

0.5 

FVR, THETA, EFP1 

72.85 

W4 

0.6 

FVR, THETA, EFP1 

65.37 

W5 

0.7 

FVR , EFP 1 , THETA 

57.83 

N6 

0.3 

FVR.EFPl, EM 5 

99.83 

N7 

0.4 

FVR.EFPl 

99.81 

N8 

0.5 

FVR.EFPl 

99.69 

N9 

0.6 

FVR.EFPl, Elf 5 

99.74 

N10 

0.7 

FVR.EFPl 

99.77 

U6 

0.3 

FVR.EFPl, WR 

99.13 

W7 

0.4 

FVR.EFPl 

98.40 

U8 

0.5 

FVR.EFPl 

98.90 

W9 

0.6 

FVR.EFPl 

99.59 

U10 

0.7 

FVR.EFPl 

99.34 
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TABLE IV- TRANSVERSE MODULUS fEE22) 


SIPPLE PIXEL 

INPUT FVR TERPE ACCEPTED 


R* 


N1 

0.3 

FVR, EFP 2 

83.50 

N2 

0.4 

FVR 

85.23 

N3 

0.5 

FVR,EFP2 

91.83 

N4 

0.6 

FVR,EFP2 

93.26 

N5 

0.7 

FVR, EFP2, THETA 

93.06 

U1 

0.3 

FVR, THETA, EFP 2 

78.36 

W2 

0.4 

FVR, THETA, EFP2 

90.73 

U3 

0.5 

FVR, THETA, EFP2 

80.15 

W4 

0.6 

FVR, THETA, EFP 2 

86.05 

W5 

0.7 

FVR, THETA, EFP2 

87.14 

N6 

0.3 

FVR, EFP 2 

87.13 

N7 

0.4 

FVR, EFP2 

86.15 

N8 

0.5 

FVR.EFP2 

90.97 

N9 

0.6 

FVR, EFP 2 

93.47 

N10 

0.7 

FVR, EFP 2 

92.05 

U6 

0.3 

FVR, EFP 2 

79.72 

W7 

0.4 

FVR, EFP 2 

70.71 

UB 

0.5 

FVR, EFP 2 

81.92 

W9 

0.6 

FVR,EFP2 

88.62 

W10 

0.7 

FVR, EFP2 

84.05 
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T£R=E V- SHEAR HXXJLOS fEC12) 


SIPPLE MODEL 


INPUT 

FUR 

TESTE flCCFPTPn 

R a 

N1 

0.3 

THETA , FVR , GPP 


N2 

0.4 

THETA , FVR , (3>p , GFP23 

97.01 

N3 

0.5 

THETA, FVR, GPP, GFP12 

98.85 

N4 

0.6 

THETA, FVR, GTF 

97.50 

N5 

0.7 

THETA, FVR, GPP, ®P1 2 

98.01 

W1 

0.3 

THETA, FVR, 9F 

98.42 

U 2 

0.4 

THETA, FVR 

94.79 

U3 

0.5 

THETA, FVR, (Tp 2 3 

94.27 

W 

0.6 

THETA, FVR 

93.71 

W5 

0.7 

THETA , FVR , GTP , GFP23 

95.62 

96.67 

N6 

0.3 

FVR, GPP 


N7 

0.4 

EVRjGPP, GFP12 

97.66 

N8 

0.5 

EVR,CTP,GFP23 

98.02 

N9 

0.6 

EVR,ap,GFP12 

96.65 

N10 

0.7 

EVR,aP,<FPi2 

97. 11 

MS 

0.3 

EVR,OP,GFP12 

98.55 

U7 

0.4 

FVR,ap,(FP12 

96.93 

MB 

0.5 

f VR,OP,GFP12 

92.45 

W9 

0.6 

fvr,qp 

95. 16 

U10 

0.7 

EVR, GPP, GFP12 

97.18 

96.90 


10 ? 


IgELE VI- POISSON’S RATIO, MAJOR (HUC19) 
SIMPLE MODEL 


INPUT 

FVR 

TERMS ACCEPTM) 

N1 

0.3 

THETA, EFP1 

N2 

0.4 

THETA, FVR 

N3 

0.5 

THETA, FVR 

N4 

0.6 

THETA, FVR, EFP1 

N5 

0.7 

THETA, FVR, EFP2 

Ml 

0.3 

THETA, EFP1 

M2 

0.4 

THETA, FVR 

M3 

0.5 

THETA 

U4 

0.6 

theta, wr 

U5 

0.7 

THETA, FVR 

N6 

0.3 

FVR 

N7 

0.4 

FVR 

N8 

0.5 

FVR 

N9 

0.6 

FVR 

N10 

0.7 

FVR 

MB 

0.3 

FVR 

M7 

0.4 

FVR, WR 

UB 

0.5 

FVR 

M9 

0.6 

FVR, GFP12,EFP2 

M10 

0.7 

FVR,EFP2 


R a 


96.39 

97.88 

96.60 

98.32 

96.62 
88.43 

84.62 

89.48 

84.05 

92.05 

97.83 

98.48 
97.77 
98.40 
99.17 
97.32 
96.45 
96.38 
98.34 
96.96 
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TAELE VII- POISSON’S RATIO, MINOR {Nnc?i ) 


SlfPLE MODEL 


INPUT 

FUR 

texie accepted 

N1 

0.3 

THETA, FUR 

N2 

0.4 

THETA, FUR, EFP1 

N3 

0.5 

THETA, FUR 

N4 

0.6 

THETA, FUR, EFP1 , EFP2 

N5 

0.7 

THETA, FVR, EFP1 

Ml 

0.3 

THETA, FUR 

W2 

0.4 

THETA, FUR, EPP2 

W3 

0.5 

THETA, FUR 

W4 

0.6 

THETA 

U5 

0.7 

THETA, FUR, EFP1 

N6 

0.3 

FUR,EFP1,EFP2 

N7 

0.4 

FVR,EFPi,EFP2 

N8 

0.5 

FUR,EFP1,EFP2 

N9 

0.6 

FVR,EFP1,EFP2 

N10 

0.7 

FUR,EFP1 ,EFP2 

M6 

0.3 

FVR, EFP1 ,GFP12 

M7 

0.4 

FVR,EFP1,EFP2 

U8 

0.5 

FVR,EFP1,EFP2 

W9 

0.6 

FVR,EFP1,EFP2 

U10 

0.7 

EFP1,FVR,EFP2 


R 2 


91.15 

94.78 

94.31 

97.18 

95.87 

90.87 
89.86 
91.93 
92.57 
94.78 

95.64 

94.90 

95.40 

93.12 

91.83 

87.73 

85.06 

84.29 

90.37 

91.42 
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INPUT 

N1 

N2 

N3 

N4 

N5 

W1 

M 2 

W3 

W4 

W5 

N6 

N7 

N8 

N9 

N10 

W6 

W7 

US 

U9 

U10 


TAg£ VIII- LONG. THERM. EXPANSION fCIElll 


SIfFLE MODEL 


FVR 

TERMS ACCEPTED 

0.3 

FVR, THETA, EFP1 

0.4 

THETA, FVR, EFPl.WR 

0.5 

FVR, THETA, EFPl.WR 

0.6 

FVR, THETA, EFP1,WR 

0.7 

THETA, FVR 

0.3 

THETA, FVR, EFP1 

0.4 

THETA, FVR 

0.5 

THETA, FVR 

0.6 

THETA, FVR, WR 

0.7 

THETA 

0.3 

FVR, EFP1 , WR 

0.4 

FVRjEFPl , WR 

0.5 

FVR,EFP1 , WR 

0.6 

fvr,efpi 

0.7 

fvr,efpi 

0.3 

FVR, EFP1 

0.4 

FVR,EFP1 

0.5 

FVR,EFP1 

0.6 

FVR.EFPl, WR 

0.7 

FVRjEFPl 


R a 


90.29 

94.46 

95.72 

95.23 

87.63 

80.53 

78.91 

84.77 

74.37 

80.50 


97.21 

96.96 

96.53 

96.60 
96.24 

91.60 
90.88 
91.55 
96.03 
94.13 
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IftHLE IK- TRAfC. TBBRM. BgfljCIQW /CTE 27 ) 


SIMPLE MODEL 


INPUT 

FVR 

terms ACCEPTED 

N1 

0.3 

FVR, THETA, WR 

N2 

0.4 

FVR, THETA, WR 

N3 

0.5 

FVR, THETA, WR 

N4 

0.6 

FVR, THETA 

N5 

0.7 

FVR, THETA 

W1 

0.3 

FVR, THETA 

U2 

0.4 

FVR, THETA, EPP1 ,WR 

U3 

0.5 

FVR, THETA 

W4 

0.6 

FVR, THETA 

U5 

0.7 

FVR, THETA 

N6 

0.3 

fvr,wr,efpi 

N7 

0.4 

FVR, WR 

N8 

0.5 

FVR, WR 

N9 

0.6 

FVR 

N10 

0.7 

FVR 

U6 

0.3 

fvr.efpi 

W7 

0.4 

FVR 

UB 

0.5 

FVR 

W9 

0.6 

FVR 

W10 

0.7 

FVR 


R s 


99.60 

99.21 

99.46 

99.69 

99.79 

95.04 

98.60 

95.19 

94.84 

97.98 

99.70 

99.53 

99.65 

99.67 

99.75 

99.15 

98.81 

98.88 

99.47 

99.22 


Ill 


TABLE X- LONG. TPEILE STRENGTH (SCXXT) 




SIfFLE I-EDEL 



INPUT 

FVR 

TERTE ACCEPTED 


R 1 

N1 

0.3 

FVR 


12.25 

N2 

0.4 

FVR, SFPT 


43.72 

N3 

0.S 

FVR 


21.68 

N4 

0.6 

FVR , SFPT, THETA 


43.68 

N5 

0.7 

FVR 


40.97 

Ml 

0.3 

FVR, SFPT 


33.37 

W2 

0.4 

SFPT, FVR 

- 

39.02 

U3 

0.5 

EFP1 , SFPT 


26.13 

W4 

0.6 

FVR.EMP 


42.27 

U5 

0.7 

FVR, SFPT 


33.55 

N6 

0.3 

FVR, SFPT 


52.12 

N7 

0.4 

FVR, SFPT, EFP1 


68.43 

N8 

0.5 

FVR, SFPT 


34.89 

N9 

0.6 

FVR, SFPT 


49.00 

N10 

0.7 

SFPT, FVR 


24.00 

M6 

0.3 

SFPT, FVR 


46.61 

W7 

0.4 

FVR 


19.33 

W8 

0.5 

FVR, SFPT 


33.13 

W9 

0.6 

FVR, SFPT 


34.40 

U10 

0.7 

SFPT, FVR 


37.65 
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TflHJ: XI- LONG. OPPRESSIVE STRENGTH (SCXKC) 

SIfFLE MODEL 

INPOT FVR TERIE ACCEPTED r 2 


HI 

0.3 

FVR 

12.25 

N2 

0.4 

FVR 

18.23 

N3 

0.5 

HONE 


N4 

0.6 

SFPC 

8.52 

N5 

0.7 

FVR 

8.08 

U1 

0.3 

WR 

8.02 

W2 

0.4 

THETA 

9.29 

U3 

0.5 

GrF,SPFC 

20.59 

W4 

0.6 

THETA 

9. 18 

U5 

0.7 

NONE 



N6 

0.3 

SFPC 

11.30 

N7 

0.4 

NONE 


N8 

0.5 

NONE 


N9 

0.6 

NONE 


N10 

0.7 

GFP12 

12.01 

US 

0.3 

FVR 

9.40 

U7 

0.4 

WR 

10.76 

U8 

0.5 

WR 

9.85 

U9 

0.6 

GFP12 

8.87 

W10 

0.7 

NONE 
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TABLE XII- TRA1SVEBSE TBBILjE STRENGTH (SCVYT) 


SIMPLE MODEL 


INPUT 

FVR 

TERTE ACCEPTED 


R s 

N1 

0.3 

FVR 


27.03 

N2 

0.4 

FVR 


32.91 

N3 

0.5 

SfFT 


8.10 

N4 

0.6 

FVR 


41.92 

N5 

0.7 

NONE 



W1 

0.3 

FVR, WR, SMPT 


26.89 

W2 

0.4 

FVR 

* 

41.43 

U3 

0.5 

FVR 


14.74 

W4 

0.6 

NONE 



U 5 

0.7 

FVR, SMPT 


31.05 


N6 

0.3 

FVR 

9.43 

N7 

0.4 

FVR 

8.19 

N8 

0.5 

FVR,EFP2 

15.58 

N9 

0.6 

NONE 


N10 

0.7 

NONE 


U6 

0.3 

FVR 

33.87 

W7 

0.4 

FVR 

13.39 

U8 

0.5 

SMPT 

8.62 

U9 

0.6 

FVR 

27.85 

U10 

0.7 

FVR 

32.77 
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SIPPLE PIXEL 


INPUT 

FVR 

terte accepted 


w* 9 

N1 

N2 

0.3 

FVR.SPFC 


R* 

33.17 

0.4 

FVR 


N3 

30.10 

0.5 

NONE 


N4 

N5 

0.6 

FVR 


38.93 

0.7 

NONE 


W1 

U2 

0.3 

FVR, WR 


28.19 

0.4 

FVR 


W3 




43.26 

0.5 

fvr,sppc 


19.57 

W4 

0.6 

NONE 


W5 

0.7 

FVR 


15.85 

N6 

0.3 

NONE 



N7 

0.4 

NONE 



N8 

0.5 

NONE 



N9 

0.6 

NONE 



N10 

0.7 

NONE 



U6 

W7 

0.3 

FVR 


28.68 

0.4 

FVR 


11.64 

MB 

0.5 

NONE 


W9 

W10 

0.6 

FVR 


31.97 

0.7 

FVR 






33.05 
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TABLE XIV- IN PLANE SHEAR STRENGTH (SOWS) 


SlfFLE MODEL 


INPUT 

FVR 

TERTE ACCEPTED 

R 1 

N1 

0.3 

FVR , THETA , GFP 1 2 

28.51 

N2 

0.4 

FVR 

8.74 

N3 

0.5 

THETA 

14.96 

N4 

0.6 

THETA, GFP 12, FVR, STPS 

31.84 

N5 

0.7 

NONE 


Ml 

0.3 

THETA ,WK, STFS, FVR 

48.16 

U2 

0.4 

FVR 

43.26 

U3 

0.5 

THETA 

8.40 

W4 

0.6 

THETA 

14.75 

M5 

0.7 

NONE 



N6 

0.3 

NONE 


N7 

0.4 

NOME 


N8 

0.5 

NONE 


N9 

0.6 

SIPS 

8.25 

N10 

0.7 

WR 

8.53 

U6 

0.3 

FVR, SIPS, GMP 

29.06 

M7 

0.4 

NONE 


U8 

0.5 

NONE 


U9 

0.6 

GFP 12, FVR 

22.20 

M10 

0.7 

sips 

17.73 
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Further regression nodels vere studied, in an attenpt to improve 
the predictive capability of the models, especially for the strengths. 
’"*** models, incorporating higher order functions and combinations of 
predictor variables used in the simple models, show some improvement 
over the simple models, proving the value of including the "interaction" 
effects of predictor variables in the regression models. In addition, 
the higher order interaction models can fit response functions over a 
wider range of fiber volume ratio, with associated improvements in the 
R 2 statistics. The data cases CONI and COM2 contain selected points 
from, the entire range of fiber volume ratios, to supply the samples for 
these runs. Furthermore, since higher order models are postulated, 

THETA ls taken to be the cosine of the average of fiber orientation 
angles. The variable fWR is a "dummy” variable, that is a function of 
other variables in the model. It is defired as 

MWR = 1 - FVR - WR 

and is intended to represent an "average" matrix volume ratio over the 

thickness of the ply. The interaction models are shown in Tables XV - 
XXVI. 

Tbe general form of the postulated models now includes higher order 

terms, so the predictor variables are tested up to the fourth power. 

Synfool ical ly , 

Y = B 0 ♦ B, (THETA) + Bj(FVR) + B 3 (WR) + B„(EFP1) + B 5 (ETP) + 
B t (MVR) + B 7 (THETA) 2 + B, (THEIA) (FVR) ♦ B, (THETA) (WR) + 

B I0 (THETA) (EFP1) + ... + B, , (THETA) *(FVR) (EFP1 ) + ... 


B 12 (THETA) M + Bj 3 (FVR ) 2 


etc. 
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The number of terms possible in a conplete fourth power polynomial 
expansion becomes unwieldy for the oases studied. Considering the 
limitation of the size of the predictor matrix in the regression package 
used (100 x 100), the terms are intuitively grouped in the hope of 
eliminating large groups at one tine. The regressions are conducted 
using "unlikely" candidates for admission into a particular ncxlel, and 
if no terms are entered, subsequent regressions are conducted without 
those terms. The justification for this approach is not a statistical 
argument, ratter an interpretation of the physical principles active in 
any chosen model. The regressions to eliminate terms are merely used as 
a check on viiat seems intuitively reasonable. 
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TABLE XV- LONGITUDINAL HXXJLUS (ECU) 


INTERACTION MODEL 


INPUT 

FVR 

TERTE ACCEPTED 


R 2 

N1 

0.3 

THETA g »FVR*EFPl 


84.50 

N2 

0.4 

TOETA g *FVK*EFPl 


92.66 

N3 

0.5 

THETA ^mFVRmEFP 1 


93.76 

N4 

0.6 

THETA g mFVRmEFP 1 


94.24 

N5 

0.7 

THETA g *FVR*EFP 1 


85.08 

Ml 

0.3 

THETA 14 mFVRmEFP 1 


63.84 

M2 

0.4 

THETA u mFVRmEFP 1 

* 

89.86 

M3 

0.5 

THETA g mFVRmEFP 1 


71.79 

M4 

0.6 

THETA 1 * mFVRmEFP 1 


64.37 

M5 

0.7 

THETA 14 mFVRmEFP 1 


55.68 

N6 

0.3 

FVRmEFP 1 , ETP 2 mIWR 


99.82 

N7 

0.4 

FVRmEFP 1 , FVR 14 


99.83 

N8 

0.5 

FVRmEFP 1 


99.72 

N9 

0.6 

FVRmEFP 1 , ETP 2 mWR , WR* 4 


99.79 

N10 

0.7 

FVRmEFP 1 , ETPmMVR 


99.79 

M6 

0.3 

FVRmEFP 1 , WR 


99.17 

M7 

0.4 

FVRmEFP 1 , MVR 2 mFVR 


98.53 

MB 

0.5 

FVRmEFP 1,WR 2 mETP 


98.99 

M9 

0.6 

FVRmEFP 1 


99.58 

M10 

0.7 

FVRmEFP 1 , ETPmMVR 


99.38 

CONI 

VARIES 

THETA ,4 mFVRmEFP 1 


96.48 

CON2 

VARIES 

FVRmEFPI, WR g 


99.92 


1X9 


INPUT FVR 


NX 

0.3 

N2 

0.4 

N3 

0.5 

N4 

0.6 

N5 

0.7 

U1 

0.3 

U 2 

0.4 

U3 

0.5 

W4 

0.6 

US 

0.7 

N6 

0.3 

N7 

0.4 

N8 

0.5 

N9 

0.6 

N10 

0.7 

US 

0.3 

W7 

0.4 

US 

0.5 

W9 

0.6 

U10 

0.7 

CONI 

VARIES 

CON2 

VARIES 


INTERACTION pedel 

TEKPE ACCEPTED 


nm*EFP2*ef. , etpz'.fvb, TOEra.Fvm.nra 

FVR*EFP2*E?F , EFP2 2 *FVR , THETA 
FVR*ETP2«ETF, THCTA 2 M£FP2,I1UR g 
FVR*EFP2#EPF , THETA 2 #EFP2 , EFP2 2 #FVR 
*** nearly singular 

THETA ^ 2 «WR , EJF 2 *EFP2 , MUR 2 
FVRkETF t THETA , EFP2 2 #FVR 
FVR*EFP2*EIF , THETA , FVR«EFP2*MVR 
THETA 2 *MUR , WR 2 #EFP2 , EJF*MUR 
*** NEARLY SINGULAR 


FVRmEFP2«EPP, EFP2 2 #FVR, E3T-MWR 

FVR*EFP2*EPF , EPP 2 FUR 
FVR*EFP2*EJF f FVR*EFP2*MVR 
FVR*EFP2wEW> , FVR 2 *QF 
FVR' 4 f EFP2«EJP 

Fra»EFP2.EH>,EFP2,EFP2’«EI*^FVR»EFF2 

FTO.ErP2«EIp t FVR«EFP2.WR,FVR.vvi 1 
FVRtfEh F2#EM? , FUR*EFP2*MUR 
FVR*EFP2*EJF , EFP2 2 «MVR 
FVR*EFP2*EPF , MUR 2 «FUR 


*** NEARLY SINGULAR 
FVR*EFP2*EMP , FVR«EFP2*MUR 


R 2 


99.19 

99.55 

98.92 

99.22 

93.26 

96.79 

93.49 

88.35 


99.22 
99.07 
98.89 
99.14 

99.23 
98.62 
98.28 
97.93 
98.44 
97.86 


99.79 
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INPUT 

N1 

N2 

N3 

N4 

N5 

Ml 

M2 

M3 

M4 

M5 

N6 

N7 

N8 

N9 

N10 

M6 

M7 

MB 

M9 

M10 

CONI 

C0N2 


TABLE XVII- IN PLANE SHEAR HDDULUS (BC12) 


INTERACTION PECEL 

FVR TERPE ACCEPTED R 2 


0.3 THETA , FVR*0*> , TffiTA ,, *FVR*G» 5 97.86 

0.4 *** NEARLY SINGULAR 

0.5 THETA 2 , FVR J #GPF I GFP12*GMP 97.75 

0.6 THETA, FVR 2 *PNR 98.01 

0.7 THETA, FVR 2 «Of > ,GFP12 98.46 

0.3 THETA* 4 , FVR ‘'♦♦Q'F, FVR 2 95.49 

0.4 *** NEARLY SINGULAR 

0.5 THETA “^FVRhPWR , THETA ‘‘wWR, WR*GPF 91.04 

0.6 THETA* 4 , FVR 3 «<JP, THETA 96.70 

0.7 *** NEARLY SINGULAR 

0.3 FVR*GPP,P1UR 2 *GFP12 97.73 

0.4 FVR*GPF,GFP12 97.97 

0.5 FVR*GrF,FVR*GFP12 96.52 

0.6 FVRwGPF,FVR 2 «<3T>12 97.10 

0.7 FVR l4 M(^,FVR»HGFP12»GPP,PWR 2 «GnP 98.90 

0.3 FVR*KW , FVR*GFP12 96.91 

0.4 FVR*GPF , FVR*GFP 1 2 92.37 

0.5 FVRK3MP,FVR«<TP12 95.08 

0.6 FVR«GMP,FVR 2 «GFP12 97.42 

0.7 FVR‘ 4 *<3 , F,FVR*<FP12*GfP,PWR 2 *CM* 96.85 

VARIES FVR 2 «WR, WR*IWR,FVR*<3 , P, THETA g *FVR»KJF 99.09 

VARIES FVR 2 <3P,WR*GPF,GFP12*<3 , F 99.54 
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INPUT 

N1 

N2 

N3 

N4 

N5 

W1 

W2 

U3 

W4 

W5 

N6 

N7 

N8 

N9 

N10 

W6 

W7 

US 

W9 

W10 


TgBLE XVIII- LONG. THERMAL EXPAfSION fCIEll) 


INTERACTION MODEL 


FVR 

terte accepted 

R 2 

0.3 

theta 2 *wr, mvr 2 , FVR*EFP1#MVB, emp 2 *efpi 

92.51 

0.4 

THETA *MVR f THETA i, ,EFP1' , ,EJF 2 *FVR,IWR 2 *EFP1 96.38 

0.5 

THETA 2 *MVE, MVR,EIF 2 *MVR, EJF 2 *EFP1 

97.26 

0.6 

FVR*EPP1 , THETA* FVR*EFP 1 , EJP 2 *MWR 

96.32 

0.7 

THETA , ETF*MVR 

90.66 

0.3 

*** NEARLY SINGULAR 


0.4 

theta m ,mvr 2 *ejf 

80.81 

0.5 

THETA EM? 2 *IWR, 

87.98 

0.6 

theta^mvr^wr 

75.20 

0.7 

theta 11 ,fvr 2 *mvr 

82.97 

0.3 

IWR 2 *ETP f EFP1 2 *EJP , FVR* 1 

99.29 

0.4 

wr 2 *emp, fvr*efpi*wr, mvr 2 *wr 

99.17 

0.5 

WR 2 *EJF, EFP1 , FUR" 

98.94 

0.6 

WR 2 *EHP , EFP1 2 *MUR 

98.94 

0.7 

wr 2 *eip ,efpi*mvr . 

99.33 

0.3 

wr 2 *etp , FVR*EFP1*MVR 

98.35 

0.4 

fWR 2 *E3V , FVR*EFP 1*MVR 

98.55 

0.5 

WR J *EIP , EFP1 2 *MUR 

98.56 

0.6 

PWR 2 *EMP , EFP 1 2 *MVR 

99.00 

0.7 

WR 2 *EFP1 , EJF 2 *MVR 

98.20 

VARIES 

THETA , fWR 3 , EFP1 2 *WR, FVR*WR*EFP1 

96.82 

VARIES 

WR 2 *EJP,FVR*EFP1*IWR | FVR M I FUR 2 *VVR. . . 

99.84 


CONI 

CON2 
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Tflm X IX- TBAjg. THERMAL Bggjgiat [ CIE22 ] 


INTERACTION MODEL 


INPUT 

FVR 

terte ACCEPTED 

R 2 

N1 

0.3 

theta 2 *wr,wr 

99.60 

N2 

0.4 

theta 2 *wr, wr 2 *fvr, fvr*efpi«em> 

99.38 

N3 

0.5 

■nffiTA 2 »WR , WR *J»FVR , FVRi 

99.48 

N4 

0.6 

theta 2 «pwr , 1WR 2 *FVR , ETP 2 »VVR, mwr 2 «vvr 

99.73 

N5 

0.1 

FVR 2 , THETA, THETA J #FVR 

99.81 

W1 

0.3 

theta 2 «wr,hvr j etpi 2 #wr 

95.16 

U 2 

0.4 

THETA 2 ***® , WR , FVR*EPP i*MVR 

98.71 

U3 

0.5 

THETA 2 *WR ( IWR 2 *FVR, EFP1 2 *WR 

95.91 

W4 

0.6 

THETA 2 ***®, THETA, THETA* 4 

95.69 

W5 

0.7 

*** NEARLY SINGULAR 


N6 

0.3 

FVRjfWR 3 

99.70 

N7 

0.4 

FVR.fWR* 4 

99.59 

NB 

0.5 

FVR.IWR 2 

99.67 

N9 

0.6 

FVR,EJP 2 »EPP1 

99.70 

N10 

0.7 

FVRjIWR* 4 

99.82 

W6 

0.3 

fvr, fvr*efpi*ejp 

99.26 

W7 

0.4 

FVR.PWR* 4 

98.97 

UB 

0.5 

FVR 

98.88 

W9 

0.6 

FVR.MVR* 4 

99.57 

W10 

0.7 

fvr,emp 2 *fvr 

99.29 

CONI 

VARIES 

theta 2 *wr 

99.32 

CON2 

VARIES 

fvr,fvr 3 ,mur 2 «eip 

99.95 
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INPUT 

HI 

N2 

N3 

N4 

N5 

W1 

U2 

U3 

W4 

U5 

N6 

N7 

N8 

N9 

N10 

U6 

W7 

MB 

U9 

W10 

CONI 

C0N2 


TABLE XX- POISSON RATIO; IlflJOB fNUC12) 


INTERACTION MODEL 


FVR 

TERMS ACCEPTED 

R* 

0.3 

*** NEARLY SINGULAR 


0.4 

THETA, EFP2*WR 

97.96 

0.3 

THETA , CFP 1 1 2*WR 

96.71 

0.6 

THETA, EFPH»rtyR 

98.17 

0.7 

THETA , FVR*EFP2 

96.48 

0.3 

*** NEARLY SINGULAR 


0.4 

THETA , THETA U *FVR*GFP1 2 

84.73 

0.5 

THETA 

89.43 

0.6 

THETA, WR*GFP 12 

B4.27 

0.7 

THETA , FVRhMVR 

92.10 

0.3 

FVR 

97.83 

0.4 

FVR 

98.48 

0.5 

FVR 

97.77 

0.6 

FVR,FVR*1WR 

98.52 

0.7 

FVR 

99.17 

0.3 

FVR 

97.32 

0.4 

FVR, WR*EFP2 

96.50 

0.5 

FVR 

96.38 

0.6 

FVR,EFP1*EFP2,GFP12*MVR 

98.41 

0.7 

FVR, FVR*EFP2 

96.97 

VARIES 

*** SINGULAR 


VARIES 

MVR, FVR*MVR,EFPl#MyR 

99.77 
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TABLE X XI- POISSON RATIO; MINOR (NUC21 ) 
INTERACTION MODEL 


INPUT 

FVR 

HHE ACCEPTED 

N1 

0.3 

THETA, THETA g *FVR*EFPl 

N2 

0.4 

THETA, FVRhEFPI 

N3 

0.5 

THETA, FVRhEFPI , EFP2*GFP12 

N4 

0.6 

THETA, THETA l **FVR*EFPl , EFP2 

N5 

0.7 

THETA , THETA g *FVR*EFPi 

U1 

0.3 

THETA, FVRmGFP 12 

M2 

0.4 

THETA , EFP2HMVR 

U3 

0.5 

THETA, FVRhEFPI 

W4 

0.6 

THETA 

W5 

0.7 

THETA , FVRhEFP 1 , THETA “hFVRhMVR , EFP2«MVR 

N6 

0.3 

FVRhEFPI , FVRhEFP2 

N7 

0.4 

FVRhEFPI , FVRHEFP2 

N8 

0.5 

FVRhEFPI , FVRHEFP2 , FVRhGFPI 2 

N9 

0.6 

FVRhEFPI , FVRhEFP 2 

N10 

0.7 

FVRhEFPI , FVRhEFP2 

U6 

0.3 

FVRhEFPI , Qrp 12 hMWR 

M7 

0.4 

FVRhEFPI , FVRHEFP2 , FVRhMVR 

ue 

0.5 

FVRhEFPI, EFP2 

U9 

0.6 

FVRhEFPI , FVRhEFP2 

W10 

0.7 

FVRhEFPI , FVRhEFP2 

CONI 

VARIES 

THETA , FRVhCFP 1 2 , EFP2 , THETA‘ , hFVRhMVR , . . . 

CON2 

VARIES 

FVRhEFPI , FVRhjwr, EFP2 , WRHGFP12 


R* 


91.69 

94.66 

95.10 

97.15 

95.82 

91.16 

89.52 
92.06 

92.53 
95.60 

95.48 

94.69 
95.52 
92.85 
91.77 

87.83 

86.48 
84.36 

89.84 
91.55 

98.70 

98.35 
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TABLE XXII- LO NGITUDINAL TENSILE STRENGTH (SCXXT) 




INTERACTION K3DEL 


INPUT 

FVR 

TERfE ACCEPTED 

R 2 

N1 

0.3 

THETA g *FVR*SFPT 

17.72 

N2 

0.4 

THETA ,, *FVR*SFPT I WR M 

47.65 

H3 

0.5 

THETA u *FVK*SFPT 

27.65 

N4 

0.6 

THETA g *FVR*SFPT, FVR 2 »EJF 

44.67 

N5 

0.7 

TRETA g *FVR«SFFT 

45.35 

Ml 

0.3 

FVK*SFPT, FVR*EFP1*HVR 

39.18 

M2 

0.4 

FVRkSFPT, fvr 2 *wr 

42.87 

M3 

0.5 

EFP 1«SFPT, EIP^IWR 

33.97 

W4 

0.6 

ErT 2 *fWR , THETA 2 *SFPT* WR 

45.09 

M5 

0.7 

THETA 2 *FVR*SFPT 

32.56 

N6 

0.3 

FVRwSFFT, FVR 2 *»ETP 

52.95 

N7 

0.4 

FVR*SFPT, FVR*EFP1 

64.41 

N8 

0.5 

FVR*SFPT, wr 2 *fvr 

39.12 

N9 

0.6 

FVR*SFPT 

47.13 

N10 

0.7 

FVR«SFPT, FVRwETF 

27.43 

U6 

0.3 

FVRwSFPT, WR 2 *EFP1 

49.71 

M 7 

0.4 

FVR*SFPT 

25.19 

U8 

0.5 

FVR*SFPT 

32.16 

W9 

0.6 

FVR*SFPT 

34.06 

W10 

0.7 

FVRwSFPT 

35.09 

CONI 

VARIES 

THETA g *FVR*SFPT, FVR«WR*WR, FVR»EFP1«MVR 

81.20 

C0N2 

VARIES 

FVRwSFPT, FVR«WR, IWR 2 «SFPT 

84.79 
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TABLE XXIII- LONGITUDINAL COW>REttQT mr 


INTERACTION MDOEL 


INPUT 

FVR 

terts ACCEPTED 

HI 

0.3 

SFPC»SfFC 

N2 

0.4 

fvr*mvr 

H3 

0.5 

NONE 

N4 

0.6 

SFPOK2V 

N5 

0.7 

FVR*GFP12 

HI 

0.3 

wr m 

U2 

0.4 

theta 14 

U3 

0.5 

etp«smpc , theta 2 *etf 

W4 

0.6 

THETA 14 

U5 

0.7 

NONE 

N6 

0.3 

SFPC , CTP 1 2*SPPC 

N7 

0.4 

NONE 

N8 

0.5 

NONE 

N9 

0.6 

NONE 

N10 

0.7 

GFP12*ETF 

W6 

0.3 

FVR*MVR 

U7 

0.4 

WR 

ue 

0.5 

WR 

U 9 

0.6 

GFP12 11 

W10 

0.7 

NONE 

CONI 

VARIES 

FVR#WR 

CON2 

VARIES 

FVR»WR, SFPC 


strength (SCXXCl 


R s 


12.53 

19.45 

9.81 

19.20 
10.49 

9.32 

23.32 

9.20 


20.04 


14.96 

11.91 

10.76 

9.85 

9.10 


46.48 

44.44 
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TABLE XXIV- TRANSVERSE TENSILE STRENGTH (SCYYT) 


INTERACTION MODEL 


INPUT 

FVR 

TERTE ACCEPTED 

N1 

0.3 

WR 2 *SfFT 

N2 

0.4 

iwr 2 *stft 

N3 

0.5 

EMPwSMPT 

N4 

0.6 

fvr 2 *iwr 

N5 

0.7 

NONE 

W 1 

0.3 

FVR 2 *WR,SfFT 

U2 

0.4 

SMPTwWR 

W3 

0.5 

mvr 2 «fvr 

W4 

0.6 

FVR*WR*EFP2 , ETP 

W5 

0.7 

FVR«SMPT 

N6 

0.3 

FVRwMVR 

N7 

0.4 

mvr 2 *fvr 

N8 

0.5 

FVR*EFP2#SI’FT 

N9 

0.6 

SMPT 2 »WR 

N10 

0.7 

FVR*EIF 

U6 

0.3 

fFR 2 *ETF 

W7 

0.4 

nm 2 *WR 

UB 

0.5 

snpr 2 *nm 

U9 

0.6 

FVR 2 *EFP2 

U10 

0.7 

mvr 2 «smpt 


VARIES THET A 11 *SI‘FT#r , WR , FVR*EFP 2*MVR , KSfFT 2 *MVR 

VARIES SfPT 2 «IWR,FVR*WR*IWR 


R 2 


31.60 

37.23 

9.61 

47.59 

25.39 

43.94 
16.32 
24. 10 
30.29 

10.47 

8.94 
13.54 

9.40 

9.13 

35.13 
19.34 
12.89 
29.27 
36.77 

73.42 

76.40 


CONI 

C0N2 
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TABLE XXV- 

TRANSVERSE OPPRESSIVE STRENGTH (SCYYCl 




INTERACTION MODEL 


INPUT 

FVR 

TERMS ACCEPTED 

R 2 

N1 

0.3 

SMPC*MVR 

33.39 

N2 

0.4 

fvr 2 *etp 

32.99 

N3 

0.5 

NONE 


N4 

0.6 

FVR 2 *MVR 

42.31 

N5 

0.7 

NONE 


W1 

0.3 

FVR 2 *WR 

26.24 

W2 

0.4 

fvr 2 *etf 

43.86 

U3 

0.5 

sw»c 2 *wr 

21.13 

W4 

0.6 

FVR*WR*EFP2 , EM 3 

25.75 

W5 

0.7 

SMPC*IWR 

18.63 


N6 

0.3 

S1FC 2 «MVR 

11.57 

N7 

0.4 

EFP2*1WR 

9.03 

N8 

0.5 

FVR*EFP2 

9.87 

N9 

0.6 

NONE 


N10 

0.7 

SM>C 2 «MVR, FVR 2 *MVR 

19.07 

U6 

0.3 

WR 2 wEPF 

32.50 

U7 

0.4 

EFP2*MVR 

14.58 

UB 

0.5 

NONE 


U9 

0.6 

fWR 2 «srpc 

32.85 

UlO 

0.7 

Mm 11 

35.79 

CONI 

VARIES 

THETA* , «S1'FC*MVR, FVR* 1 

76.43 

C0N2 

VARIES 

IWR‘ I ,FVR 2 «1WR I MVR 2 «SMPC 

75.59 
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T QELE XXVI- IN PLA N E SHEAR STENGTH fSTXVg} 


INTERACTION MODEL 


INPUT 

FVR 

TERfE ACCEPTED 

N1 

0.3 

FVRwGFPl 2*K3T , THETA‘S 

N2 

0.4 

FVR*GFP12#EIF 

N3 

0.5 

THETA 

N4 

0.6 

1HETA ,, *GFPi2,SIPSwMVR 

N5 

0.7 

NONE 

Ml 

0.3 

THETA, FVRwVVRwEJF, THETA a wSrFS, FVRwMVR 

U2 

0.4 

THETA g wFVR,THETA g wCFP12 

W3 

0.5 

THETA g 

W4 

0.6 

THETA, FVR* WR 

W5 

0.7 

THETA **»FVK 

N6 

0.3 

NONE 

N7 

0.4 

SIPSwMVR 

N8 

0.5 

NONE 

N9 

0.6 

SlPSjSMPS 1 * 

N10 

0.7 

fvr 2 *pwr 

US 

0.3 

SMPSwMVR.GMP’* 

W7 

0.4 

FVR*GFP12*IWR 

H8 

0.5 

NONE 

U9 

0.6 

FVRwGFPl 2*IWR 

U10 

0.7 

SfPS 

CONI 

VARIES 

THETA g «FVR , FVR 2 wSIPS 

C0N2 

VARIES 

FVRwWR,MVR , *,FVR a wVVR 


R J 


27.64 

13.51 

14.97 

30.84 

52.20 

26.58 

12.89 

22.33 

10.72 


11.24 

16.14 

11.40 

28.58 

8.28 

19.20 

17.73 

36.74 
61.46 


CHAPTER IV 


DISCUSSION 


A. Owruiew 

U* minerical simulations conducted show that certain assumptions 
about the statistical distribution of local nonuniformities in fiber 
composites lead directly to gtiantif iable uariations in material 
properties. The advantages inherent in the stochastic characterization 
are numerous. The deuelopment of quality control and reliabilty 
measures for composites is crucial to their acceptance in aircraft 
designs. The rediction in needed experimental data achieoahle throujh 
jtmlicious simulation of the wide uariety of available composite material 
systems could significantly lo»r the costs of material selection and 
acceptance testing. In the results of this stilly, the confidence 
intervals calculated can be interpreted as the predict of an 


emperimental program, specifically designs as an analog of the physical 
processes vfcich occur in real materials. 
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®* Histograms and Distributions 

Data cases 1, 2, and 3 demonstrate the differences between a 
deterministic base case and random cases with narrow and wide dispersion 
of input data about the base case. 

In Fig. 30, it is apparent that the deterministic case 1 vali* of 
15750 ksi. for longitudinal nodulus falls near the nean of the case 2 
data. However, the case 3 sanple appears to have a nean slightly lower 
(approximately 15000 ksi.). It should be noted that the size of the 

interval over vhich the sanple occurs is noticeably larger in the widely 
distributed case 3 run. 

Transverse modulus, (Fig. 31) demonstrates a higter nean value for 
the wide distribution than for the narrow, which is greater than the 
deterministic value of 1065 ksi. reported in Table II. The increased 
transverse modulus is related to the added stiffness available in fibers 
with high misalignnent relative to longitudinal direction. 

Shear modulus, (Fig. 32) is measurably changed by nonuniformities. 
Hie deterministic value of 516 ksi is exceeded by the case 2 valie of 
approxinately 620 ksi, which is further exceeded by the case 3 value 

near 900 ksi. Fiber misalignent has a significant effect in shear 
vrodulus variation. 

Poisson’s ratios (Fig. 33, 34) show similar trends in location of 
sanple means and relative dispersion of the sanple for the data studied. 

Poisson’s ratios generally increase with fiber misalignent and volume 
fraction changes. 
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IT* coefficients of thermal expansion (Figs. 35, 36) for the sample 
studied reflect the longitudinal contraction of graphite fibers «*>en 
heated. The longitudinal coefficient of thermal expansion for 
AS-graphite fiber is -0.550 x 1 <T‘/ F, while the transverse coefficient 
is 0.560 x 10 / F. The offset orientation of crystal lattice planes 

in graphite fibers can explain this behavior. These values, the fiber 
misalignment, and fiber volume ratio near 0.5 all contribute to the 
occurrence of a negative longitudinal coefficient of thermal expansion 
for the conposite. At higher fiber volume ratios, the values calculated 
vould be less than in the present case, because of the contolling fiber 
behavior for high fiber volune ratio. 

The longitudinal strengths (Fig. 38, 39) are significantly reduced 
V *' en nonunii r ° rmi ties are present. The deterministic case 1 value of 203 
ksi. for tensile strength is compared to a mean near 160 ksi for case 2 
and a mean near 130 for case 3. In compression, the deterministic value 
of 165 ksi. compares to means near 100 ksi. and 80 ksi. for the narrow 
and wide distributions, respectively. The failure mode in compression 
varies in the random sanples . 

Transverse strengths (Fig. 40, 41) show sensitivity to the 
variations assumed. Misalignments, volume fraction nonuniformities, and 
constituent strength variations all contibute to reduction in the 
strength val«s. Sub-ply shear failures occur, which undermine the 
already low transverse composite strengths. 

In plane shear strength (Fig. 42) values decline from 10.01 ksi. 
for case 1 to a mean near 8.0 ksi. for case 2. However, case 3 shows a 
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^ — 8 0 *'“• 14 -B— ‘ha, , te added 5heir , tlTOgth 

d- to fiber misalign i, balanced by the r«l„oed stre ^ th , q 

variable fiber volume fraction. 


C. Confidence Curves 

Tte effects of carious shape parameters of fiber strength are sho™, 
“ Figs. 43 and 11. The higher wibull distribution shape parameter of 
» prodt^es a narrow distribution of fiber strength values. The 
composite that has few weaker fibers is expected to be. stronger, and 
Fig. 13 demonstrates this few lonitudinal tensile strength. However, 
compressit* failure (Fig. 41) is a more complex phe-xmenon. I„ the 
region of low fiber volume ratio, the -rule of .natures' failure 
criteria for a subply can control the failure mode, at higher fiber 
uolume ratio, however, compressive failure can be initiated by 
delamination, or by a shear failure in a sub-ply. The mixt.ee of 
failure modes in compressive failure is not .ell understood, but can 
explain the seeming inconsistency of the intersection of the curves in 
Fig. 44. at a fiber volume of 8.7, the weakest fibers (a . 18) are in 

the strongest composite, when strength is normalized with resp^t to 
fiber conpressive strength. 

nte effects of various shape parameters for matrix stre^ths are 
staled in Figs. 45, 46, and 47. Transverse tensile a^ compressive 
Strengths show expected redtotions for lower matrix strengths. I„- P l m 
shear strength show, lover dispers.cn a, a large fiber voUm. of 8.7, 
and also declines in general for higher fiber volume. 
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The flber misalignment effects are studied in Figs. 48-57. 
Longitudinal modulus (Fig. 48) shows narrow intervals and slight 
reductions for greater misalignment. Transverse nodulus (Fig. 49) and 
in plane shear nodulus (Fig. 50) are enhanced by fiber misalignment. 
Longitudinal tensile and conpressiw strengths are degraded by 
misalignment (Figs. 51, 52). Transverse tensile aid conpressive 
strengths are enhanced (Figs. 53, 54). In-plane stear strength shows 
total separation of confidence intervals between curves with different 
degrees of misalignnent . Poisson’s ratios (Figs. 56, 57) increase for 
high fiber misalignnent values. 

The fiber stiffness effects (Figs. 58-67) are very snail for the 
distribution paraneters studied. 

D. Examination of Regression Models 

The regression models for thernoelastic properties demonstrate 
resonably high predictive capability in the simple models assuned. 
Iferginal improvement s are achieved in expanding the models to include 
higher order interaction terms. Further improvement is gained by using 
sanple data from the wide range of volume percent valiss. Tfc? higher 
multiple correlation coefficients of these models may be die to the 
additional information available in the sample size of 100 that was 
used. The nearly singular predictor matrices which occur in the higher 
order models indicate that terms must by selectively removed to 
eliminate linearity between assuned predictor terms. The regression 
results support the use of the simple models for ttermoelastic 
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properties, because ieprouenents in prediction capability i» the higher 
order „odel S for the Mre data are ,n»ll. 

Strengths are not modeled wll by the slnple the intepaoUon 
»deu. The predictors chosen are a«rage preppies, vfereas «» 
strengths are based on the seahest points in the xaterial. Eden the 
unidirectional cases (NS-N 1 ®, OS-MO, present data that the interaction 
"Odels ha„e considers, e difficulty in accoeodating. Sosevhat gre„. r 
predictiue uaioe is gained by using the expanded data for strength nodel 
prediction. Osin, fourth order algebraic functions, ualues of the 
luiltiple correlation coefficient square approach 857. for longitixiinal 
tensiie strength. The other strength, generally base poorer result,. 


CHAPTER V 


CONCLUSIONS 


A tractable, constituent based, probabilistic analysis procedure 
for fiber conposites has been developed using the ICAN program as a 
basis. Within the limitations of the mechanics of .material model, 
properties and strengths of a variety of composite material 
conf igurat ions can be sinxilated. 

m>is stirty quantifies the thernoelastic and strength properties of 
a graphite/epoxy ply subject to assured uncertainties for fiber 
misalignnent t constituent volune fractions, and constituent properties. 
He* results show seueral advantages of probabilistic characterization of 
this xaterial. These include the identification of unforseen variations 
in conposite material properties, and the nechanical effects of local 
nonuniformities. The relative inporta™, of the various fabrication and 
material variables on conposite properties is identified, and the 
resulting behavior quantified. 

advantages of a probabilistic formulation of conposite naterial 
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properties ooer a deterministic one are numerous. Comparison of the 
results of this Study with test data could reveal some sources of 
previously unaccounted scatter in the data. Expected val« ranges could 
be predicted for experimental results. Since the situations provide 
data that is analagous to experimental data at lover cost, laboratory 
classification, material selection, and acceptance testing of composites 
can be guided by the infornetion nude available by these methods. 

Although the method presented provides results for only the basic 
Ply, emtension of the simulation to incite lamination .angle variations 
in a general layup is feasible. Since finite element material property 
cards are generated, structural analysis of components with randomly 
varied properties defined at a number of points in the body can supply a 

"=re realistic description of the random nature of stcpta-al response 
due to material inhonogeneity. 


•Hie stochastic formation of material properties is generally 
recognized as one requirenent of failure theories for materials, 
ftlthoujh the failure criteria in the rodels used in this study are 
conservative, progressive failure of fiber conposites could be node led 
by incorporating load redistribution and material property recalculation 
in the vicinity of failed material. 
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This appendix outlines the theories and equations in the ICflN 
program that are used in this project. In tte first section on 
conposite micromechanics, the elastic and thernal properties of a 
conposite ply are defined with respect to its principal naterial axes 
The next section, dew ted to laminate theory, contains the 
transformations and sunnations of ply properties used to arrive at 
laminate properties. The last section contains a brief discussion of 
the failure criteria. 


1. Conposite micronechanics 

The theory for calculation of the properties of a unidirectional 
fiber conposite ply based on the properties, volume fractions, and 
orientation of its constitwnts is knovn as conposite micromechanics. 

In this section, the subscripts f, m, v, and / represent fiber, natrix, 
void, and laminate, respectively. The synbolic notation and the 
equations used are sunnarized below: 

Volune fractions: 


k + k ♦ k =1 

F m v 


Longitudinal fbdulus: 


E /ll * k f E fll + k m E m 


Transverse Mbdulus : 
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E /22 “ E /33 = 


1 ~ ^7 (1 ~ E /E„„) 
* ' m £22 7 


Shear Itoduli: 


i - v^r (i - g /g i 
* m £12' 


1 ~ v£7 (1 - G /G ) 
* m £23 7 


Poisson’s Ratios: 


^12 - W /13 “ + k /( v /i 2 ~ •») 


U /23 = k *> 'Vo-1 + k 2u - — £ 

f /23 n m E /n *722 


^ oe ^^'^ c ^ en 't s of thernal expansion 

a r 1 1 + k T [a E /E ) 
Pll m Ll mm £11 ‘ 

a /ll = 

1 + k «< E . /E m - >) 




“/ 22 - 


-E— r 


1 * V. E m 


E rn + l „< E n - I * a t22 k t 


a 33 = a /22 
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2. Laminate Tteory 

■Hus section describes the methods viiich are used to calculate the 
elastic properties of laminates from the properties, orientation, and 
distribution of individual laminae. The elastic properties are then 
used to predict the response of the laminate to external loads. The 
nethods used to predict stresses in the laminae under application of 
external loads are also described. Failure loads can be predicted by 
using these methods; as described in a following section. 

a. Generalized Hooke’s Law 

The stresses acting at a point in a solid can be represented by the 
stresses acting on the planes noriml to the coordinate directions, or 
equivalently, on the surfaces of an infinitesimal cube as showi in Fig. 
B-l. The stresses (o.^) on each face are resolved into three 
conponents: one nornel stress and two shearing stresses. The first 
subscript refers to the direction nornel to the plane in vrfiich the 
stress acts and the second subscript to the direction in viiich the 
stress acts. The stress conponents shown on the faces of the cube are 
taken as positive and can be taken as the forces (per unit area) exerted 
by the naterial outside the cube upon the material inside. A stress 
conponent is positive if it acts in the positive direction on a positive 
face of the cube. Thus normal tensile stresses are positive, and nornal 
conpressive stresses are negative. Nine stress conponents mist be used 
to define the state of stress at a point, nanely ^ ^ 

°31» a l2* a 32’ a !3’ a 2L‘ ’ n,ere **■ nine corresponding strain 


1S8 


components, following the sane subscript convention. 

For bodies in vfcich each strain conponent is a linear function of 
all six stress components, the generalized Hooke’s Law can be expressed 

a ij - E ijkl f kl 

'* ere E ijkl is a fourth order tensor of elastic constants. For nine 

stress components and nine strain components, there mist be 81 elastic 

constants defining E^^. Certain reductions in tf» nuriber of 

independent constants for an anisotropic body are due to symretry 

properties of the tensor E^^. By considering nonent equilibrium about 

the center of the cube, it can be shown that at any point a = a 

23 32’ 

°31 = a i3’ and °12 = a 21 ‘ E ijkl is syimetric with respect to the 

first tvc indices. Second, because the strains are syimetric (that is, 

f ij = f ji ^ ’ E ijkl nust ^ syn^Etric with respect to the second two 
indices. This reduces the number of elastic constants to 36. Furtter 
reduction to the final 21 elastic constants for a gereral anisotropic 
material is accomplished by assuming the existerce of a strain energy 
density function, such that 

U = U{s. .) 

with the property 

au 

af »j “ ° 1J 

From the generalized Hooke’s Law, 

au 

df. . ~ ijkl f kl 

differentiation with respect to yields 
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— MSLl - E 

de kl i J ijkl 

Since the order of partial differentiation is immaterial, 


a 

fJ5L] 

a 

f au J 

de ki 

la*. .] 

ij J 

= ds. . 

krJ 


and the subscripts can be interchanged to yield 



so that 


E ijkl “ ^lij 

the first pair of subscripts in can be interchanged with tte 

second pair without any change in the val*s. The nunber of elastic 
constants is thus reduced to 21. 

b. Lamina Constitutive Relation 

Several siirplif ications to the generalized Ifooke's Law can be Trade 

for the special case of a thin orthotropic material, vhich approximates 

a unidirectional fiber composite lamina. By considering the invariance 

of elastic properties under coordinate transformation for planes of 

symmetry, the tensor E..^ can be rediced to the following nine 
constants: 


E im 

E 1122 

E 1133 

E 1 122 

E 

2222 

E 2233 

E 1133 

E 

2233 

E 

3333 


0 

0 

1 0 

It is now convenient to make the following notation changes: 
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a ll 

= a i 

°22 

“ °2 

a 33 

= a 3 

°23 

= T 23 

°13 

" T 13 

a l2 

" T 12 


*11 = *1 


22 


= o„ 


33 * 3 
2 ®23 = y 23 = *4 
2S 13 = V 13 = *5 
2e !2 = y !2 = *6 


The generalized form of Hooke's Law can now be written 


6 

L 

</=l 


a i ~ ^ C . .5 . f or 1 , j r I C 

1 .. J,/ J A * • • • f b 


matrix C ijF is known as the stiffness retrix, andV are the 

J 


a l 


*2 


a 3 

II 

T 23 


CJ 


KJ 



nponents. 

In 

retrix 

form 

Hooke' 

s Law j 

C 11 

C 12 

C 13 

0 

0 

0 


f i 

C 12 

C 22 

C 23 

0 

0 

0 


*2 

C 13 

C 23 

C 33 

0 

0 

0 


*3 

0 

0 

0 

C 44 

0 

0 


y 23 

0 

0 

0 

0' 

c„ 

0 


y 





35 



31 

0 

0 

0 

0 

0 < 

- 


y 






"66 


12 


material. For laminae that are assured sufficiently thin, the through 

the thickness stresses are zero. Thus * 3 = <* 4 = ° 5 = 0, for plane 

stress. It is apparent that s =5=0 

4 5 

T*e stress strain relations for a thin unidirectional lamina are 


written 
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*11 *12 • 


'I 

^12 ®22 0 


*2 



1 

» ® ^66 


T V 12 j 


l 

using the tensorial strain y V ^2 instead of the engineering strain 
The Q terms are known as reduced stiffnesses, i.e. 


Ei 

0 - r - 

V 11 11 1 - V i 2 V 2 i 

v 12 E 2 v 2 ]Ej 

12 12 1 - U] 2 u 2 | 1 • v 1 2 u 2 1 

E 2 

*22 = C 22 = 1 - o l2 u 21 

*66 = T ^ C ll ” C 12^ = G 12 

v * iere Ej, E 2i v i2 , w 2 i, and Gj 2 are the ply elastic constants, neasured 
with respect to the natural iraterial system. It ney be noted that only 
four of these constants are independent. 

The stress- strain relation above shows that there is no coupling 
betveen tensile and shear strains, as long as the applied stresses are 
coincident with the principal neterial directions. However, coupling 
appears vften a lamina is tested at arbitrary angles with respect to the 
principal nsterial directions. The general form of tte stress-strain 
relation for any angular orientation of a lamina is considered next, 
c. Stiffness matrix transformations 

A lamina is loaded along a coordinate system x-y oriented at sone 
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angle # with respect to the principal material directions as show in 
Fig. B-2. Since stress and strain are seoond order tensors, they are 
transforned by 


and 


s, 

i 


°1 


a 

X 

°2 

= m 

a 

y 

T 1 2 J 


T 

*y } 

1 

1 

s 

X 

2 

*12 

- [T] 

£ 

y 

i 

tv i 


xy ' 


vhere [T] is the transformation matrix for plane stress and plane strain 

transforned by clockwise rotation about the (3,z) axes, given by 

sin 2 # 2 sin# cos# 

cos 2 # -2 sin# cos# 


[T] = 


cos 2 # 

sin 2 # 


~sin# cos# sin# cos# cos 2 # - sin 2 # 


Inversion and substitution yields 


a 

X 


5 

X 

a 

y 


s 

y 



i 

T 

xy J 


r * 

xy' 


vhich is the stress strain relation for a lamina referred to arbitrary 
axes. For sinplicity, the notation [ Q ] is introduced 

[Q] = [T] _1 [Q][T] 

viiere [Q] is called the transforned reduced stiffness iratrix. 

Using the approach outlined above, it is possible to obtain 
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expressions for the elastic properties referred to the x-y coordinate 
system. 

Elastic properties of laminates 

A number of assumptions are trade in laminate theory to obtain 
theoretical predictions. Ttese are: 



1 . 

2 . 

3. 


the lamina are perfectly bonded and do not 
each other 


slip relative to 


?° nd betveen the laminae is inf initesiirally thin 
the laminate has the properties of a thin steet 


■mese assumptions allow the laminate to be treated as a thin 
elastic plate. The classical hypothesis of Kirchhoff is applied to 
derive the strain distribution throughout the plate under external 
Forces, fecause the laminate is composed of laminae oriented in 
different directions with respect to each other, the stress-strain 
equation for each layer (k) is defined as 

1 f *11 *12 *16 1 f * x 

°Y = *12 *22 *26 E y 

T *y I *16 *26 *66 ik l T V xy 

"Huis for a given strain distribution, the stress in each layer can be 

determined. The strain at any point in a laminate undergoing 

deformation must be related to the displacements and curvatures of its 

midplane. The discussion vfcich follows assumes that the laminate is 

thin. Kirchhoff plate theory is used in this formulation. 

TTie deformation of an arbitrary section of a laminate is show in 
Fig. B-3. It is assumed that lines straight and perpendicular to the 
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" UdPla " e tof0re deEorra “°„ rein so otter defornmtion. This is 

•quivelent to negloo, in, transoms. shewing derations. Conning 
Fig. B-4(b) with Fi,. in *, ich t „ t „ 

r«Bin perpendicular after defection, it is seen that the upper and 

loser surfaces of the plate rust no, shift their relative positions. I, 

is ohoious that the resistance of a thin plate to such defection is 

large, rech larger than its resistance to defections perpendicular to 
the midplane. 

It is assuned that the point B a, the midplane undergoes 
displaceirents u 0 , u 0 , and „ 0 along the », y, and z a,es, respectively. 
The displacenent u in the „ direction of a point C located on the nocl 
flBCD at a distance z from the midplane is given by 

u = u 0 - za 

vtere a is the slope of the midplane in the x direction, 


a = 


dw 0 

dx 


The last tvo equations can be used to obtain the displacenent u of an 
arbitrary point at a distance s from tlx? midplane as 

dw 0 


Similarly, 


u = u 0 - 2 


v = v 0 - z 


dx 


dw 0 

~dy 


Since the strains noml to the midplane are neglected (plane 

strain), the displacement w at any point is taken equal to the 

displacenent w 0 at the midplane. The strains in terms of displacenent u 
and v are 
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du 

du 0 

d 2 w® 

x “ Ox 

Ox 

" 2 ~0^ 

dv 

dv 0 

d 2 w 0 

y = dy c 

W 


du 

dv 

du 0 i 


xy 


dv 0 d 2 w 0 

dy + ax = dy~ + dx 2z dxdy 


In terms of midp lane strains and plate curvatures, tte strains in a 
laminate vary linearly through the thickness, 
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viiere midplane strains are given by 


f° i 


du 0 

x 


dx 



dv 0 
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dy 

y° 


du 0 dv 0 

xy J 


; dy dx 


and the plate curvatures by 
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a 2 u 0 


dx 2 
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dxdy 


The stresses in any (k) lamina can be obtained by substituting the 
previous equation into the stress strain equation 

Q, 
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e. Laminate Stiffness 1-fetrix 

Classical laminate theory provides a method of oaloulatii^ the 
resultant forces and moments per unit length acting on tte laminate by 
integrating the stresses acting in each lamina through tte thickress (h) 
of the laminate. Resultant forces are obtained by 


N = 
x 


h/2 


-h/2 


a dz 
x 


N 

y 


f h/2 


-h/2 


a dz 
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N 

xy 


h/2 


-h/2 


T 

xy 


dz 


The TiDnent resultants are obtained by integration through the thickness 
of the corresponding mo Rents of stresses about the midplane: 


M 


x 


M 



T h/2 

J “h/2 

T h/2 

J -h/2 

r h/2 

J -h/2 


o' z dz 

X 


a z dz 

y 


T 

xy 


z dz 


The units of N, N, N are force per unit length and H . H . M are 
x y xy x y xy 

mment per unit length. The sign conventions are shown in Fig. B-5. 

Using the resultant force and mment relations, a system is defined 
that is statically equivalent to the laminate stress system, but applied 
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at the midp lane. Thus, the external loading has been reduced to a 
system that does not contain the laminate thickness or z coordinate 
explicitly. 

For a laminate consisting of n laminae (Fig. B-6), the resultant 
force-monent system acting at the midplane can be obtained by adding 
integrals representing the contribution of each layer by 




Using the expressions for the stresses in the k-th lamina derived 


earlier, and noting that the mi dp lane strains and plate curvatures are 
constant not only within the lamina, but for all laminae, it is apparent 
that they can be taken outside the integral sign. Tlx? stiffness netrix 
[Q] is constant within a lamina so it also can be taken outside tte 
integration to give 
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Three new net rices, A. .. B . 
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i are def ined , Were 


A. . 

ij 


I < V 


k th k " ^t-i) 


k=l 


B. . = 
ij 


D. . = 

ij 


1 

~ 2 1 V* **4 - ti) 

k=l 
1 " 

— 2 < Vk I’i - *i-i> 

k=l 


These new natrices, A, B, and D, simplify the resultant force and moment 
relations, and are know as the extensional, coupling, and bending 

stiffness matrices, respectively. The total plate constitutive equation 
is then 



It nay be recalled that in an orthotropic lamina with arbitrary 
orientation the shear stress is coupled with the normal strain and the 
normal stresses are coupled with the shear strain. In general, a 
resultant shearing force on a laminated plate produces midplane normal 
strains in addition to the expected shearing strain. Similarly, a 
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resultant normal force will induce shear strains in addition to midplane 
normal strains. 

The nonzero coupling matrix B in the plate constitutive equation 
explains the coupling between bending and extension of the laminated 
plate. Thus, normal and shear forces at the midplane induce not only 
midplane de for rat ions, (and hence, midplane strains) but also twistirx? 
and bending, producing plate curvatures. Similarly, resultant bending 
and twisting moments induce midplane strains. 

f . Lamina stresses and strains 

The aim of the analysis of a laminated conposite is to determine 
the stresses and strains in each of the laminae forming the laminate. 
These stresses and strains are used with failure criteria to predict the 
loads for failure initiation for a laminate. The failure criteria are 
discussed in the section devoted specifically to that purpose. 

The strains in a lamina caused by external loading are a function 
of laminate midplane strains and plate curvatures, as previously 
discussed. Once the lamina strains are known, lamina stresses can be 
found using the lamina stress-strain law. Thus, the starting point for 
calculating lamina stresses is the determination of laminate midplane 
strains and plate curvatures in terms of the applied loading. The plate 
const itutiw equation given previously can be inverted to give the 
midplane strains and plate curvatures explicitly in terms of the 
resultant external forces and ncments. The result of the inversion 
process is 



I * j l C f D f 

J ’ er ' ! B’, and D’ are simplified forms of the inversion process 
results, and are functions of the A, B, and D matrices of the original 
form of the plate constitutive equation. 

It is no» apparent that math these equations, an analysis of a 
laminate subjected to external forces and moments can be conduced: 

1. calculate midplane strains and plate curvatures 

,.:i - 

calculate lamina stresses in global ( K -y) system 
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calculate lamina stresses in natural 
transverse to fiber) svstem. 


(longitudinal arri 
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The strain variations in a lamina are calculated in an analagous 
■numer. The stress-strain variation is compared wth the alloMhle 
stresses and strains in each lamina. Thus the load at .Inch failure is 
initiated in one of the lamina can be calculated, as long as a strength 
criteria exists in terms of the lamina natural Mis system. The 
formulation of lamina failure criteria is discussed in the next section. 
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3. Strength Tteories 

It is .ssurol that the strength of a laminate „, t t „ 

^ Str '” 9tM ° f «" i " di '" d “* 1— •. • .1*1. failure criteria 

consists of evalisting the lamina strengths in their principai snteria, 

directions subject to induced stresses or strains at the boundaries of 
the lamina. In this context, it is assumed that the lamina and its 
constittmnts behaue in a linear elastic nanner to failure. The strength 
analysis described here assures that the behavior of each lamina in an 
arbitrary laminate is the sane as the behavior observed in the natural 
axis system when the lamina is part of any o, te r lami«te urrfer sane 
resses or strains. In other words, it is assisted that the strength 
criteria for a lamina i„ plane stress is valid for any orientation of 
the lamina in a laminate. In the ICON program, the lamina strengths are 
calculated using the expressions given below. 

Longitudinal tension 

S I11T * S fT <*f * WW 
Longitixiinal conpression: 

The longitudinal conpressioe strength .mist be conputed on the basis 
of three different criteria: 

a. rule of mixtures 

s /nc - s «: ( k f ♦ k „vw 

b- delamination 


3 /llC 


* (13 S /12 * S vl 
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c. fiber microbuck ling 


F- O 
2 iff 


‘/lie 1 - VI " 


'Transverse tension 


S ,22I ' S„ T (F*CT/tO»l) 


Transverse compression 


S = S„ / DENOtt 
j22C mC 


Transverse shear 
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v*iere Fj and F 2 are given by 


n 

r i = J “4k" 


f 2 = 1 - 


4k 


■ 

^ m 


The variable EEHOM is introduced for convenience: 


Donn = [1 - ^7(1 - e b /e„ 2 )3 rrrfFT-ry-T V,,r - !)■ 


vtere V is given by 


F 1 ■ W 1 - -V - VWi 


T - 


F 1 ‘ 
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The variable FACT is used to correlate the strengths of HTE and Kevlar 

oonposites with the experinentally observed values. Since neitter 
of these fibers is used in this work, FACT takes the value unity. 
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Fig. B.3- Bending geometry in the x-z plane. 
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